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Abstract 

We investigate Voronoi-like tessellations of bipartite quadrangulations on surfaces of arbi- 
trary genus, by using a natural generalization of a bijection of Marcus and Schaeffer allowing 
one to encode such structures by labeled maps with a fixed number of faces. We investi- 
gate the scaling limits of the latter. Applications include asymptotic enumeration results 
for quadrangulations, and typical metric properties of randomly sampled quadrangulations. 
In particular, we show that scaling limits of these random quadrangulations are such that 
almost every pair of points are linked by a unique geodesic. 

Nous examinons les proprietes de mosaiques de type Voronoi sur des quadrangulations 
biparties de genre arbitraire. Ceci est rendu possible par unc generalisation naturcUe d'une bi- 
jection de Marcus ct Schaeffer, permcttant de dccrirc ccs mosaiques par des cartes ctiquetees 
avec un nombrc fixe de faces, dont nous dctcrminons les limites d'echelle. Parmi les applica- 
tions de ces rcsultats, figurcnt Ic comptagc asymptotiquc des quadrangulations, ainsi que des 
proprietes mctriqucs typiques de quadrangulations choisies au hasard. En particulicr. nous 
montrons que les limites d'echelles de quadrangulations aleatoires sont telles que presque 
toute paire de points est liee par un unique chemin geodesique. 

Keywords. Random maps, scaling limits, random snakes, asymptotic enumeration, geodesies. 
MSG. 60C05 ; 05C30 ; 60F05 

1 Introduction 
1.1 Motivation 

A map is a graph embedded into a compact orientable surface without boundary, yielding a 
cell decomposition of the surface, and considered up to orientation-preserving homeomorphisms. 
Random maps are considered in the physics literature on quantum gravity as discretized versions 
of random surfaces [2]. This approach allows one to perform computations of certain integrals 
with respect to an ill-defined measure on surfaces, by approximating them by finite sums over 
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maps. From a mathematical perspective, this leads to the stimulating problem of existence of a 
measure on compact surfaces arising as the scaling limit of, say, uniform random triangulations 
of the sphere (or the compact orientable surface of genus g) with a large number of faces. 

Important progress has been made in this direction in the recent years, thanks to bijective 
approaches initiated in Schaeffer's thesis [12]. They allowed Chassaing & Schaeffer [11] to obtain 
the scaling limits for the radius and profile of a uniform rooted planar quadrangulation with n 
faces, considered as a metric space by endowing its vertices with the usual graph distance. In 
particular, they showed that typical graph distances between vertices of such random quadran- 
gulations are of order n^^^. Generalizations of this result have been obtained for much more 
general families of random maps by Marckert, Weill and the author in [60\ [5^ 1^5} I36j . relying 
on generalizations of Schaeffer's bijection by Bouttier, Di Francesco &; Guitter [7]. 

An important step has been accomplished by Le Gall [26], who showed that scaling limits of 
these random quadr angulations, considered in the Gromov-Hausdorff sense [22], must be home- 
omorphic to a metric quotient of the so-called Brownian Continuum Random Tree of Aldous [1] . 
This random quotient was introduced under the name of Brownian map by Marckert & Mokka- 
dem [32], who proved convergence of random quadrangulations towards this limit in a sense 
yet different from convergence in distribution in the Gromov-Hausdorff topology. Subsequently, 
Le Gall & Paulin [28] inferred that the Brownian map is homeomorphic to the two-dimensional 
sphere. At the present stage, it is however not known whether the scaling limit is uniquely deter- 
mined as the Brownian map, which would lead to a satisfactory answer to the above mentioned 
problem. 

A natural idea is to characterize the limit through its "finite-dimensional marginal distri- 
butions". From the point of view of metric measure spaces, this would mean to characterize 
the laws of mutual distances between an arbitrary number of randomly sampled points. This, 
however, seems hard to obtain. The approach of the present work, which deals with combina- 
torial, probabilistic and geometric aspects of maps, was initially motivated by another notion 
of finite-dimensional marginals, which takes into account the planar structure of the graphs 
that are considered, and roughly consists in studying the Voronoi tessellation on the map with 
sources taken at random. 

Our central combinatorial result is a natural generalization of Schaeffer's bijective construc- 
tion, to be introduced in Section [2j We do not restrict ourselves to the planar the most 
natural framework is to consider maps of arbitrary (but fixed) genus g, so our bijection is really 
a generalized version of the Marcus-Schaeffer bijection in arbitrary genus [SStllOj. Let us explain 
briefly the idea of our construction. Recall that the Marcus-Schaeffer bijection encodes a genus-^f 
bipartite quadrangulation with a distinguished vertex by a labeled map of genus g with one face, 
in such a way that the labels keep track of graph distances to the distinguished vertex in the 
quadrangulation. Our construction encodes a genus-(j( bipartite quadrangulation with k distin- 
guished vertices by a labeled map of genus g with k faces (Theorem U] and Corollary [1]) . In some 
sense, the faces of the encoding map correspond to the Voronoi tiles in the quadrangulation, 
with sources at the k distinguished vertices, and the labels allow to keep track of graph distance 
of a vertex to the source of the Voronoi tile it belongs to. An important fact is that there is 
some flexibility in the construction, allowing one to consider a generalized Voronoi tessellation 
in which the distances to the k sources are measured after addition of a delay depending on the 
source. 

The resulting set of labeled maps with a fixed number of faces is combinatorially much 
simpler than that of bipartite quadrangulations, and we are able to study their scaling limits. 
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This is done in Section U] (Theorem [5] and Propositions [3ll3|), after showing how to encode labeled 
maps by appropriate processes in Section [3l Sections [5] and [7] give three applications of these 
scaling limit results for labeled maps, to combinatorial and geometric properties of bipartite 
quadrangulations 

The first application is a new derivation of known asymptotic enumeration results (Theorem 
d]), to which the short Section [5] is devoted. These are initially due to Bender &: Canfield [1], 
who obtained the asymptotic number of rooted maps of genus g with n edge^ by recursive 
decomposition methods, of a very different nature from our bijective study. These results have 
also been obtained in the recent work of Chapuy, Marcus & Schaeffer |10] . who completed the 
exact enumeration of maps of genus g that was initiated in . The starting point of our study 
(using the Marcus- Schaeffer bijection, which corresponds to the case A: = 1 of our bijection) is 
the same as in [10], but the rest of our approach is different, the arguments of Sect. [J] being of 
a more probabilistic nature. 

Our other results concern the metric structure of randomly sampled bipartite quadrangu- 
lations. Let us indicate briefly what we mean by "random" in this paper. Most the articles 
on the topic have focused on scaling limits as n — > oo of random families of maps conditioned 
to have n faces or n vertices. We prefer to randomize the number of faces, by using natural 
cr-finite measures on bipartite quadrangulations of fixed genus g, called the Boltzmann measures 
Qg (the term being inspired from [3], see also [301 ED), and which are obtained by assigning 
an appropriate weight to the faces. In the Physics terminology, these measures correspond to 
the so-called grand- canonical measures, while the measures with fixed number of faces are the 
microcanonical measures. These measures are used to define natural "Boltzmann-Gibbs" dis- 
tributions on quadrangulations, depending on an inverse temperature parameter (3 that allows 
one to tune the average size of the quadrangulation and take scaling limits. 

The existence of these scaling limits, considered with respect to the Gromov-Hausdorff- 
Prokhorov topology on metric measure spaces, and stated in Theorem [21 is then obtained from 
the study of Section HI This generalizes to arbitrary genera the fact that planar quadrangulations 
admit scaling limits for the Gromov-Hausdorff topology, as shown in Le Gall [2^. The theory 
of metrics on weighted metric spaces that is needed here is developed in Section [6l 

The main result of the present paper. Theorem [31 gives qualitative information on the metric 
structure of the scaling limits of random bipartite quadrangulations. We show that these scaling 
limits are geodesic weighted metric spaces, in which two typical points are linked by a unique 
geodesic. The idea is to use our bijection for quadrangulations with k = 2 marked vertices, in 
order to show that, ii x,y are typical points in the limiting random metric space (X, d) arising 
as a scaling limit of random bipartite quadrangulations, the intersection of a geodesic between 
X and y and a geometric locus of the form 

{z e X : d(x, z) - d{y, z) = D} , D G [-d{x, y),d{x, y)] , 

is a.s. reduced to a single point. 

Theorem [3| is strongly reminiscent of the combinatorial results of the recent paper by Bouttier 
& Guitter [S]. They show that in the planar case g = 0, two typical vertices in a random 
quadrangulation with n faces cannot be linked by two geodesies which are "very far" in the 
scale n^/^, as n — > oo. Our result confirms that geodesies between typical points become unique 

^This is equivalent to our result, since such maps are in one-to-one correspondence with bipartite quadrangu- 
lations of same genus with n faces [331 Proposition 1]. 
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in the scaling limit, answering one of the questions raised in Section 5 of [8], though in our slightly 
different setting of Boltzmann-distributed maps. We also mention the recent work by Le Gall 
[27] on similar topics as the present paper and [8], and providing results on the exceptional 
geodesies as well. 

In this work, the cardinality of a set A is denoted by |^|. 

1.2 Embedded graphs and maps 

Let us introduce some formalism for embedded graphs and maps, which we partly borrow from 

m- 

Graphs on surfaces. Let 5 be a compact connected orientable two-dimensional surface with- 
out boundary. It is well-known that such surfaces are characterized, up to homeomorphism, by 
an integer g > 0, called the genus of S, and that the topology of S is that of the connected sum 
Sg of g tori. The surface of genus is the two-dimensional sphere. 

A half-edge, or oriented edge e in 5 is a continuous path c : [0, 1] — > S, considered up to 
continuous increasing reparametrization, which is either injective on [0, 1], or is injective on [0, 1) 
and satisfies c(0) = c(l). In the latter case, e is called a loop. The order on [0,1] induces a 
natural orientation of the half-edge e, and its start and end points are defined as e~ = c(0) and 
e"*" = c(l), a definition that depends only on the half-edge e and not on the function c. Similarly, 
the image of e is defined as Im(e) = c([0, 1]), and the interior of e is c((0, 1)). The reversal e of 
the half-edge e is defined as the reparametrization class of c(l — •), so that e"*" = e~ ,e~ = e"*". 
An edge is a pair of the form e = {e, e} where e is a half-edge. 

An embedded graph in S, or simply a graph on S, is a pair {V,E), where 

• E is a non-empty, finite set of edges such that two distinct edges intersect, if at all, only 
at some of their endpoints 

• y is the set of vertices, i.e. of end-points of elements of E. 

We let E be the set of half-edges corresponding to E. The set E has even cardinality and 
comes with the involution e i— > e. An orientation of the edges is a choice of one half-edge inside 
each edge to form a set -£^1/2 C E with same cardinality as E. 

A graph on S determines faces, which are the connected components of the complement of 
V U UeeB 5'. Let F be the set of faces of the graph. 

Maps. We say that the embedded graph (V, E) is a map if all faces are simply connected (hence 
determining a cell complex structure on S). A map is necessarily a connected graph by |29t 
Lemma 1.5]. 

The degree of a face in a map is the number of edges it is incident to, where it is understood 
that edges that are incident twice to the same face should be counted twice. By Euler's formula, 

\V\-\B\ + \F\=x{g), (1) 

where xid) = 2 — is the Euler characteristic of S. 

A map is rooted if it comes with a distinguished half-edge e* G E. 

We declare two (rooted) maps m, m' on surfaces S, S' to be isomorphic if there exists a 
homeomorphism S ^ S' that preserves the orientation, and sends the vertices and edges of m 
to those of m' (and maps the root of m to that of m'). Note that the number and degrees of 
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vertices, faces, edges and genus are preserved under such transformations. We will systematically 
identify maps that are isomorphic, making the set of maps a countable set. 

We will usually denote (equivalence classes of) maps by bold lowercase letters m, q, . . .. The 
standard notation for the root of a rooted map will be e*, and it should be clear according to 
the context to which map e* refers. Also, when dealing with maps, it will often be understood 
that we are working with a particular graph belonging to the class of that map. This allows to 
consider the set of (half-)edges, vertices and faces of a (isomorphism class of a) map m, denoted 
by £^(m), E(m), V{m), F{m), or simply E, E, V, F when there is no ambiguity. 

A quadrangulation is a map whose faces all have degree 4. A map is called bipartite if its 
vertices can be colored in black or white, in such a way that two neighboring vertices are assigned 
different colors. We let Qg be the set of rooted bipartite quadrangulation of genus g (planar 
quadrangulations are all bipartite, but it is not the case in higher genus). For every q G Qg we 
have 4|F(q)| = 2|E(q)| = |F(q)| and |y(q)| = |F(q)| +x(5) by ©. 

Distances in maps. For a map m, and two vertices x, y, we say that the half-edges ei . . . , 6^ 
form a chain from x to y if = x,e^ = y and ef = e^j^ for 1 < i < n — 1. We define the 
graph distance dm{-c,y) as the minimal n such that there exists a chain of n half-edges from x 
to y, and a chain from x to y with length dm{x,y) is called geodesic. 

1.3 Spaces of (weighted) metric spaces 

The main motivation of this paper is to study the metric structure of large random maps, 
endowed with their graph distance. This needs to introduce some concepts of metric geometry, 
essentially originating from Gromov's work [22], which have developed a growing interest in the 
probabilist community, starting from works of Evans and coauthors. See |16 t ll7 [ [ ^ ll5 t l26j. All 
results that are needed here will be presented in Section El 

A weighted metric space is a triple {X, d, n) where {X, d) is metric space and fi a Borel 
probability distribution on {X,d). Two metric spaces are isometric if there exists a bijective 
isometry from X onto X' . The isometry class of the space {X,d) is denoted by [X, d]. Two 
weighted metric spaces {X, d, fj.) and {X' , d' , fj,') are said to be isometry-equivalent if there exists 
a bijective isometry (p : X ^ X' such that = fi' . The isometry-equivalence class of {X,d,fi) 
is denoted by [X, d, fi] . 

We let M be the set of isometry classes of compact metric spaces, and the set of isometry- 
equivalence classes of compact weighted metric spaces. As we will see in Section El it is possible 
to endow M„ with a distance cIghP; the Gromov-Hausdorfi'-Prokhorov distance, that makes 
(Mw,dGHp) a Polish space, see Theorem El and Proposition El The topology on will be 
implicitly that inherited from cIghp- If a > 0, and X= [X,d] eM (resp. X = [X,d,fi] G M^), 
we let aX = [X,ad] (resp. aX = [X,ad,fj]). 

Also, a metric space {X, d) is called a geodesic metric space if for which every pair of points 
x,y G X, there exists a geodesic path between x and y, i.e. an isometry 'jxy : [0,d{x,y)] — > X, 
with 7a;j/(0) = x,'^xy{d{x^y)) = y. In particular, geodesic metric spaces are arcwise connected. 
For compact metric spaces, being a geodesic space is equivalent to being a path metric space, 
i.e. an arcwise connected metric space (X, d) such that the distance d[x, y) is achieved as the 
minimum over all continuous paths c : [0, 1] ^ X with c(0) = x,c(l) = y, of the length of c, 



1 INTRODUCTION 



6 



defined as 

fn-l 



length(c) = sup 1^ (i(c(ti),c(ti+i)) : n > 1, = to < ti < . . . < tn = l| • 

We let PM be the set of isometry classes of compact path (or geodesic) metric spaces. 
1.4 Main results 

Our core combinatorial result (Theorem H]) is a generalization of the Marcus-Schaeffer bijection, 
encoding fe-pointed bipartite quadrangulations (with delays between sources) with the help of 
labeled maps with k faces. The central statement in this paper is Theorem [5] (and its exceptions, 
Propositions [3] and Hj), which allows us to take scaling limits of this simpler class of maps. We do 
not reproduce the exact statements of these theorems in this paragraph, as important notational 
background is needed. 

We give several applications of these results to various (combinatorial and probabilistic) 
aspects of maps. The first is a counting result, giving the following asymptotic behavior for the 
cardinality of the set of bipartite quadrangulations of genus g with n faces, consistently with 

mm- 

Theorem 1 For any g >0, it holds that 

|Q"i ~ C„ 12" n-5>^(^')/^ 

for some constant Cg £ (0, oo) defined at \21\) . 

In the case 5 = 0, this result is a trivial consequence of the exact formula (see [11] ) 

2 

Qo ~ — r~n 2 Cat„ , 
n + 1 

where Cat^ = ;^^(^^) is the n-th Catalan number, yielding the value Co = 27r~^/^. 

From this statement, we see that the generating function for bipartite quadrangulations 
counted with respect to the number of faces has radius of convergence 12~^, and it is natural 
to introduce the critical Boltzmann measure Qg on assigning mass 12~l'^(*')l to the element 
q. Since |l^(q)| = |-F'(q)| + x(ff); the function /? i— > Qg(e~^^'i) is analytic in /3 G (0,oo), where 

= |y(q)| is the "volume" of q. We define the Boltzmann- Gibbs probability distribution on Qg 
with inverse temperature P by 



,,,_ Q,(e-^^..dq) 
^« - Q,(e-/^^.) ' 



and 

QgiV^e-^^'^) ■ 

The reason for the difference between the cases 5 = and 5 > 1 is technical and will be explained 
in due course. 
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With every q G Qg we associate the metric space Xq = [y(q),(iq] and the weighted space 
Xq = [y(q),(iq,;Liq], whcrc 

is the uniform distribution on the vertices of q. 

Theorem 2 For every g > and P > 0, the probability measures 

QW-^[{{a'y'X-,a-'V^)e-}), a>l (2) 

form a relatively compact family of probability distributions on M„ x (0,cx3), endowed with the 
topology of weak convergence. 

Theorem 3 A limiting point S^g^^ of ([^, which is a probability measure on M„ x (0,oo), is 
supported on spaces {[X,d, ^],V) such that 

1. [X,dL\ is an element o/PM, 

2. is a diffuse measure of full support, 

3. for /i (8> /i almost every x,y £ X, there exists a unique geodesic joining x and y. 
Comments. 

1. Introducing the extra parameter P and dividing it by a parameter a ^ oo is referred to in 
Physics as "approaching the critical point" . One possible interpretation is to see it as a device to 
compute various Laplace transforms of "continuum" limits for discrete models (see also Remark 
3.). It is also reminiscent of the basic method of singularity analysis consisting in approaching 
a singular point, and this is exactly what we do in Sect. [5] as a special case. On a very basic 
level, one should just note that letting a — > oo has the effect of making the expectation under 
Qlf^°'^ of quantities such as go to infinity, so that we are intuitively taking a limit as the size 
of the map gets large. Alternatively, we could also have considered conditioned laws of the form 
Qg{- I {K~^ < a^^Vq < K}) for some fixed K > 1, and let a — > oo. 

2. Considering uniform measures on bipartite quadrangulations with n faces amounts to con- 
ditioning the measure Qif^ on {|i^(q)| = n}, regardless of (3, and a similar result as Theorems 
[2] and [3] could thus be considered as a conditioned version of these statements. We believe that 
such generalizations can be obtained at the price of technical complications. Note that in the 
case g = 0, a. result similar to Theorem [2] appears in Le Gall [26] in such a conditioned set- 
ting, though the convergence does not really take /.iq into account (see however \2()\ Remark a., 
p.646]). 

3. As mentioned in the beginning of the article, a satisfactory statement would be to get rid of 
the extraction in Theorem O by showing that there is actually a unique limiting point to the 
set of measures ([2]). A further guess would be that this unique measure could be written in the 
form =^g(e-^^d(X, V))/=^g(e"^^) for some cj-finite measure Yg on x (0,cx)). 

4. Last, a natural question is to ask whether our results are universal. That is, is it still true 
for, say, Boltzmann-Gibbs distributions on maps of genus g with faces of arbitrary valences 
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and weight wi^. on a face of degree fc? The natural generalization of Qg would be the measure 
assigning weight Y\ /GF(m) ^deg^ / on the rooted map m, for a particular "critical" tuning of the 
weights w. In this vein, the results of \30\ [3H HSj [36] should be generalizable to our context. 

Acknowledgments. This work was initiated while visiting the Centre for Mathematical Studies 
in Cambridge (UK), supported by a London Mathematical Society funding, where the question 
about uniqueness of geodesies arose in motivating discussions with James Norris. Thanks also 
to Guillaume Chapuy and Gilles Schaeffer for interesting discussions around [inj- We finally 
thank the two anonymous referees and the Editorial Board for very useful comments, which 
were greatly appreciated. 

2 Delayed Voronoi tessellations of bipartite quadrangulations 
2.1 Fat graphs 

As mentioned before, maps are of combinatorial nature, and there exists a purely combinatorial 
description of maps which will be useful to our purpose. A nice introduction to the different 
definitions of maps can be found in |24[ Chapter 1]. 

A map m can be encoded in a triple of permutation^. Start from a particular representation 
of the map m on the surface S. Around each vertex, the outgoing half-edges are cyclically 
ordered in a non-ambiguous way, by considering the first intersection points of the half-edges 
with a small circle around the vertex, oriented counterclockwise. These cyclic orders associated 
with distinct vertices involve distinct elements of E, and thus determine a permutation cJm of 
E, whose cycles are naturally associated with the vertex set V. We also let am(e) = e, and 

= (^m'^m- ^he triple (cTm, am, ^m) is a fatgraph structure on E. We will usually drop the 
subscripts m when the situation is unambiguous. 

To understand what is, the reader might check on examples that it corresponds to the 
following. Visiting the boundary of each face of the graph in counterclockwise order, determines 
unambiguously a cyclic ordering of the half-edges belonging to this oriented boundary. The cycles 
of (fm are exactly the so-defined cyclic orders, which are in a one-to-one natural correspondence 
with faces of the graph. See Figured] for an illustration. One deduces that the group generated 
by (frrn ctm, '/5m) acts transitively on E, amounting to the fact that maps are connected. 

Conversely, every triple {a,a,ip) of permutations on a finite set E such that ipaa = 1, a 
is an involution without fixed points, and the permutation group generated by (cr, a, (p) acts 
transitively on E, can be represented in the above way as the triple associated to a uniquely 
determined map of some genus g, determined by Euler's formula. 

A fancy intuitive image of a fatgraph is to consider it as a road network with cars driving 
on the left lane. An edge and its reversal constitute the two lanes of the road, while the vertices 
are roundabouts connecting the roads. The cycles of a correspond to the outgoing roads from 
the roundabouts, the cycles of (p correspond to the roads taken by a driver always taking first 
left at each roundabout, and the cycles of a correspond to a U-turn. 

In a fatgraph, each half-edge e belongs to exactly one cycle of a and ip, and thus determines 
unambiguously a vertex and face of the map, called the vertex and face that are incident to e: 
in terms of graphs, they correspond to the origin vertex of e and the face lying to the left of 



^Note that in [24], the permutations act to the right, while we prefer to make them act to the left, which results 
in differences like the formula ipaa = 1 below instead of aaip = 1, referenced in Remark 1.3.19 in this book. 
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Figure 1: A portion of a map, emphasizing a face of degree 7, and an illustration of the fatgraph 
structure. The action of the permutations a, a, (p are shown on a particular half-edge e. 



e, respectively. By slight abuse of notation, for / a face of a map, we will write e £ / if the 
half-edge e is incident to /, i.e. belongs to the cycle corresponding to / in the fatgraph structure. 

If 61,62 are two distinct half-edges incident to the same face, we can break the cycle of ip 
containing 61,62 into two linearly ordered sequences of edges which are denoted by [61,62] and 
[62,61]. Formally: 



so that [e,e] = {e} with this definition. We also define (61,62] = [61,62] \ {61}, and so on. 

Another notion that will be useful is that of corner. In a map, the corner incident to the 
half-edge e is the intersection of a small neighborhood around the vertex incident to e, and the 
angular region comprised between e and a{e) (to be perfectly accurate, we should define the 
corner as the germ of such sets). In the sequel, we will often speak indistinctly of a half-edge or 
its incident corner. 

2.2 Quadrangulations with sources and delays 

We first consider unrooted maps, and will re-inject the root afterwards. For g > 0,k > 1, we let 
Q° ^ be the set of 3-tuples of the form (q, x, D) where 

(Ql) q is a bipartite (unrooted) quadrangulation of genus g, 

(Q2) X = (xi, . . . ,Xk) are elements of V{q), 

(Q3) D = [di,...,dfc] is an element of the set P(q, x) of delay vectors, i.e. elements of Z'^ 
defined modulo addition of a common integer, such that 



[61, 62] = {p'{ei) : i > 0, 62 ^ {^{ei), < j < i - I}} 



for every 



1 < i / j < 



(3) 
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and such that 

dq{xi,Xj) + di - dj £ , l<i,j<k. (4) 

When k = 1 this definition is void, and D = [0] will usually be omitted from the notation. Note 
that a delay vector such as in (Q3) exists if and only if miuj^j (iq(xj, Xj) > 2. Indeed, in that 
case we can take a bicoloration of V^(q) in black and white, and notice that two vertices are at 
even distance if and only if they are of the same color. Setting = if Xi is white, and dj = 1 if 
Xi is black, yields a vector (di, . . . , dfc) with the wanted property. Conversely, if d{xi,Xj) = 1 for 
some i ^ j, then obviously we should have di ^ dj and \di — dj\ < 1 which is impossible. This 
implies some restrictions on the quadrangulations involved in Q° ^. — for instance, they should 
have at least k distinct vertices. 

With any element of Q° ^, we can associate a partition of the set of edges of q that is very 
much reminiscent of Voronoi tessellations in a Euclidean space. We fix arbitrarily a represen- 
tative of a delay vector {di, . . . ,dk), and define li{x) = df^{x,Xi) + di for x S V{(l). Then for 
1 < i < k, and for x and y adjacent in q, 

\h{x)-h{y)\ = l, (5) 

since this quantity can be at most one and vertices are at even distance of Xi if and only if they 
have the same color in a bipartite coloration of the vertices. Thus, defining 

l{x) = mm{di + df^{x,Xi)) , 

i 

we deduce that |l(x) — l(y)| < 1 for every adjacent x,y. 

On the other hand, by (Q3) we have d{xi, Xi)+di — di € 2Z, which shows that li{xi) — lj{xi) € 
2Z for every By dSJ, this property propagates from xi to all other vertices, and since the 
graph is connected it holds that Ij — Ij takes its values in 2Z for every Consequently, we 
cannot have li{x) = lj{y) with adjacent x,y for any and this implies that l{x) ^ l{y) for 
adjacent x,y. Therefore 

— l(y)| = 1 , for all adjacent x, y . 

This allows us to fix an orientation E1/2 C E{q) by enforcing 1(6"*") = l(e~) — 1 for all 
e G -£^1/2- The function 1 decreases along chains of edges of -Ei/21 &iid plainly, the only vertices 
without edges pointing outwards are xi, . . . ,Xk- Thus, any chain of edges of £'1/2 can be extended 
into a maximal chain ending at some element of {xi, . . . , Xk}- Also, any such chain with initial 
vertex x is geodesic, and the endpoint of a maximal extension must be a Xi with i such that 
l(x) = li{x). To see this, recall that each of the functions Ij makes a step in {—1, 1} when passing 
from a vertex to one of its neighbors, by Q. Moreover, we just saw that the function 1 decreases 
by 1 when going along one step in the chain. So if the chain terminates at Xj, it must hold that 
l{y) = li{y) for every vertex y along the chain. Since l(x) = (iq(xi,x) + di and l(xj) = di, the 
chain has length (iq(xj,x) and is a geodesic. 

We color each edge of £'1/2 with one in k colors in the following way. For e S -E1/2, consider 
the maximal path starting with e and turning as much to the left as possible, that is, after 
arriving at a vertex, it takes the first outgoing edge it encounters in clockwise order around the 
vertex. If this path ends at vertex Xj, assign the color i to the edge e. 

Definition 1 The set £'j(q, x, Z?) of edges colored i is called the D-delayed Voronoi tile with 
sources x and center Xj. 



2 DELAYED VORONOI TESSELLATIONS 



11 



The D-delayed Voronoi tiles form a partition of -Bi/2- The reason for the name is that when 
D = 0, the set £'j(q, x, 0) contains all the edges pointing from a vertex that is strictly closer to Xi 
than any other source, which is the usual property of Voronoi tiles in Euclidean space. A general 
delay vector D corresponds to the fact that distances from different sources are measured with 
a relative advantage of di — dj of source xj on source Xi. Note that it might be that no edge 
pointing from a vertex x for which dq{x, Xi) = miuj (iq(x, Xj) lies in -Ej(q, x, 0). In fact, it holds 
that if 

Vi{q_, X, D) = {xi] U {e- : e G Ei(q, x, D)} , 1 < i < A; , 
then these sets cover F(q) and we have 

|x G V{q) : li{x) = l(x) < m_inlj(x)| C Vi{q,x,D) C {x G V{q) : li{x) = l(x)} . 

It is not difficult to see that the first inclusion is always strict as soon as /c > 2, because one can 
find some x G Vi(q, x, D) and some j ^ i such that x G Vj{q,x, D), and therefore it holds that 
lj(x) = Ij(x) = l(x). The second inclusion can be strict, although this is harder to see at first: on 
FigureEl the only vertex x with label 3 is at distance 3 of xi and 2 of X2, and di = 0,d2 = 1, so 
that Ii(x) = l2(x) = l(x), however, x does not belong to Vi{q, x, D), because the (only) leftmost 
geodesic started from x ends at X2- 

2.3 The bijection 

An embedded graph is said to be labeled if its vertices are assigned integer values 1 = (I(x),x G 
V), in such a way that |l(x) — l(y)| < 1 for every adjacent x,y. Labeled maps are defined 
accordingly. We let LM° ^ be the set of pairs of the form (m, [1]) such that 

(LMl) m is a map of genus g with k faces indexed as /i(m), . . . , /fc(in) 

(LM2) 1 = (l(x),x G V{m)) is a labeling of m, and [1] is its class up to an additive constant in 
Z. 

Our goal in the remaining of this section is to construct and study a bijective mapping 
k between Q° ^ and LM° ;., and its inverse mapping We first make some preliminary 

remarks. 

1. The construction starts from a particular embedding of an element of Q° ^, and builds an 
embedding of the image labeled map, by deletion of some of the edges and vertices of the initial 
quadr angulation, and addition of new edges. 

In particular, if (q, x, Z?) and (m, [1]) correspond by the construction, the set of vertices of 
m will be identified with the set V{q) \ {xi, . . . ,Xfc} of vertices of q distinct from the sources. 
This will be important in the sequel. This is also the reason why we will give the same name 1 
to labeling functions on a quadrangulation q and its image map m, as the labeling function on 
m is indeed the restriction to the vertices of m of the labeling function 1 on q introduced above. 

2. The mapping \I'° . is a natural analog of a bijection of Schaeffer between pointed planar maps 
and well-labeled trees [HIE]. More precisely, in the case {g,k) = (0,1), the set LMq ^ is con- 
stituted of labeled (unrooted) planar trees, which themselves are in one-to-one correspondence 
with well-labeled trees, for which all labels are positive and at least one vertex has label 1, by 
taking any labeling 1 and adding — mini + 1. A generalization of Schaeffer's bijection to genus 
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Figure 2: An example of a 2-pointed quadrangulation of genus 0, with delays di = 0,d2 = 1 
(left) and the associated labeled map with 2 faces by \I'q 2 

g > 1 is discussed in Marcus & Schaeffer |33j . It is recovered as the special case /c = 1 of our 
construction, where the set LM° is constituted of labeled maps of genus g with only one face, 
called g -trees in [33] . 

3. We can understand our generalized bijection as running simultaneously the Marcus-Schaeffer 
construction at k distinct competing vertices. Informally, let water flow at unit speed from the 
sources xi, . . . in such a way that the water starts diffusing from Xi at time di, and takes 
unit time to go through an edge. When water currents emanating from different edges meet at 
a vertex (whenever the water initially comes from the same source of from different sources), 
they can go on flowing into unvisited edges only by respecting the rules of a roundabout, i.e. 
edges that can be attained by turning around the vertex counterclockwise and not crossing any 
other current. The process ends when the water cannot flow anymore, and the tile Ei{c[,x,D) 
is the set of reversed half-edges that are visited by the water flowing from source Xj. The parity 
condition on D implies that the water flows emanating from different sources never meet in the 
middle of an edge, but always at vertices. 

We now give the rigorous construction, and refer to Figure [2] for an example. 
2.4 Construction of the bijection 

We start from an element (q, x, D) G Q° ^, and take some particular embedding of q. Choosing 
arbitrarily (di, . . . ,dk) G D, we label the vertices by the function \ : x G V{q) 1-^ minj((ij + 
(iq(x, Xi)). By ([5]), the labels of vertices encountered when exploring a face of q, starting from the 
vertex with lowest label and turning counterclockwise, must be either of the form I, l + l, 1+2, l+l 
or / + 1, Z + 1. Following [33] , we call faces simple and confluent accordingly. We then perform 
the construction of |33] on this labeled object, that is, we add extra edges to the map m, one 
inside each face: 
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/ + 1 



1 + 2 





l + l 



l + l 



Figure 3: Adding new edges (thick lines) inside simple and confluent faces of a graphical repre- 
sentation of an element of Q° ^ 



• if the face is simple, then the added edge links the vertex with highest label to its successor 
in clockwise order around the face, hence splitting the face into a face with degree 4, and 
a face of degree 2 

• if the face is confluent, then the added edge links the two vertices with highest label, hence 
splitting the face into two triangular faces. 

In this way, we have deflned a new map (graph) q' whose faces have degree 2, 3 or 4, and with 
labeled vertices, whose set is still V{q). See Figure El 

The newly added edges, excluding the ones of the original map, define a graph m on the 
surface on which q is drawn. It is easily obtained from the construction that the vertex set of 
this new graph is exactly y(q)\{xi, . . . , Xk}- Indeed, the vertices xi, . . . ,Xk have neighbors with 
strictly larger labels, and thus have the least labels in any face they are incident to, and cannot 
be incident to a newly added edge. Conversely, if x is a vertex not in {xi, . . . ,Xfc}, with label 
l(x) = I, we let ei, . . . , be the outgoing edge cycle associated with x in the fatgraph structure 
of q. Then there exists i such that l(e^) = 1 — 1. It is then easy that the new edge added in the 
face of q incident to Cj (i.e. located to the right of Cj) must be incident to x, regardless of the 
face being simple or confluent. 

Lemma 1 The embedded graph m is a map. Each face of m contains exactly one element of 
{xi, . . . , Xk}, and is indexed accordingly as /i(m), . . . , /^(m). 

Once this is proved, and recalling that vertices of q are vertices of m as well, we use 1 = 
(l(x),x G ^(q) \ {xi,... ,Xfc}), as a labeling function on V{m.). It is now straightforward to 
check that ^'°_^(q, x, L>) := (m, [1]) satisfies properties (LM1,LM2). 

Proof of Lemma [H We proceed in a similar way as in [33]. Namely, we consider the dual 
edges of q' that do not cross the newly added edges, and give them an orientation according 
to Figure HJ i.e. in such a way that the vertex located to the right of the oriented dual edge 
has strictly smaller label than the one located to the left. In doing so, we create an oriented 
spanning subgraph of the dual graph of q', call E'^^^ t^ie set of its edges. 

By Figure HI for every face / of q', there is a single element of Ey^ that goes out of it. 
Therefore, there is a unique, infinite chain of E'^,^, say (e^,m > 1), starting from /. This 
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Figure 4: Drawing and orienting the dual edges not crossing the newly added edges inside each 
of the possible types of faces of q' 



infinite chain eventually cycles. By inspection of Figure HI if /j is the label of the vertex of 
q located to the right of ej^, then it must hold that > Im+i for every m > 1. But since 
(e^,m > 1) eventually cycles, the sequence {lm,m > 1) is eventually constant, equal to / say. 
Another look at the first and third situations of Figure [3] shows that the only possibility is that 
the cycle (i.e. the sequence e^,m > N for n large enough) surrounds a single vertex of q. This 
vertex must be an element of the set of sources {xi, . . . , x^}, since these are the only vertices of 
q' that are incident only to edges of q, and hence that can be surrounded by a cycle of E'^^^- 

This implies that the (oriented) graph E'^^^ has k distinct connected components, and hence 
m has k faces. Now starting from the union of faces and edges of q' that are incident to Xi, 
which is a topological disk because Xi is incident exactly once to each face, we can add the faces 
of q' leading to Xi in the dual path of edges of E°^2 without ever making a cycle, so that the 
faces of m are simply connected, hence m is a map. □ 



2.5 Converse construction 

In the sequel, the label 1(e) of a half-edge (or the associated corner) e will be a notational 
shorthand for the label of its origin l(e~). 

The idea of the converse construction is simply to perform the Marcus-Schaeffer inverse 
construction separately inside each face of m. Specifically, given (m, [1]), we add an extra 
vertex Xi inside the i-th. face of m. We choose an arbitrary labeling 1 G [1] and extend it to 
V{m)U{xi,...,Xk} by 

l(xj) = min 1(e) — 1 . 

e£fi{m):e-=x 

We let (cj, a, ip) be the permutation 3-tuple associated with the fatgraph structure on m. 

We define the successor s{e) of an edge e E fiii^) by s(e) = Xi if 1(e) = l{xi) + 1, and 
otherwise, 

Ke)=¥,™(^)(e), 

where 

m(e) = min{m > 1 : l(v3"(e)) = 1(e) - 1} . 
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We then draw an edge between the corners incident to each edges e and s(e) for each e G E(m), 
in such a way that these edges do not cross each other. This is possible, because by the 
constraint on the labels, we have |l((/9'"(e)) — l((/?'""'"^(e))| < 1, which entails that if e has successor 
s(e) ^ {xi, . . . , Xk} and e' G [e, s(e)], then we have l(e') > 1(e) and thus s(e') must be in (e', s(e)]. 
If s{e) = Xi then e has minimal label inside the face it is incident to, and the labeling constraint 
prevents e from being in any set of the form [e' , s{e')]. 

When the new edges have been added in this proper non-crossing way, we delete the edges 
of m and keep only the newly added ones. 

Lemma 2 The following properties hold. 

1. The graph thus obtained is a map q, which is a bipartite quadrangulation. 

2. The labels G ^(q)) are equal to = mini<j<fc((iq(xi, x) + di), where di = \{xi). 

3. Setting D = [di, . . . , dk], the triple (q, x, D) is an element of Q° ^. 

4- The set of oriented edges of q from e to s{e), for e £ /i(m), is £'.j(q, x, D). 

Proof. We first check property 1. Let e be an edge of m with label 1(e) = /. Its successor 
s(e) has label / — 1. If ip{e) has also label /, then its successor s{ip{e)) = s(e), so that the edges 
between e, (/^(e) and their common successor form a triangle with the base edge e. This is the 
same in the face incident to the reversal e, so that the deletion of {e,e} after adding the new 
edges yields a 4-valent face. 

If (/j(e) 7^ I, then either ^{e) has label / + 1, or ip{e) = s(e) by definition of the successor. In 
the first case, notice that we also have that s{s{ip{e))) = s(e), because s{ip{e)) must have label 
I and lie in [e,s(e)]. Thus, the newly added edges between (the corners associated with) e and 
s(e), (p{e) and s{ip{e)), and s((^(e)) and s(e), must form a 4-valent face with the edge e. Now, 
in the face e is incident to, observe that the successor of e must be a~^{e) = (fie), as they are 
successive vertices with labels / + 1 and I. Thus, the newly added edges around the edges e,e 
form a 4-valent face. It remains to treat the case where s(e) = ip{e), but then we are in the 
situation just discussed in the face incident to e. 

Therefore, any (non-oriented) edge of m is surrounded by a 4-valent face in the new graph, 
and these are the only possible faces, as any new edge is drawn from some corner of m. Thus, 
the new graph is a map on the same surface as m. It has only square faces, and coloring vertices 
according to the parity of their labels yields a bipartite coloration. 

To check property 2., let x be a vertex of m incident to the face /j(m), with label I, and 
let e G /j(m) be an edge pointing out from x. Then the (half)-edges drawn between the 
corners e, s(e), s^(e), . . . , form a chain with length l(x) — l(xj), and by definition, it turns 
to the left as much as possible among the half-edges of q along which 1 decreases. Moreover, 
any path using only the edges in q must be at least as long, because these edges always link 
vertices with labels that differ by 1 exactly. Hence, this path is a geodesic and one deduces 
dq{xi,v) + di = I. For the same reason, we must have dq{xj,v) + dj > I for every j G {1, . . . , k}, 
so that I = mmj{dq{v, Xj) + dj). This gives, in passing, that the edge of q from e to s{e) is an 
element of £'j(q, x, D), yielding the property 4. 

We now check that D = [di, . . . ,dk] is a delay vector, which will entail property 3. First, 
we have seen that for any vertex x, the length of a shortest path from Xi to Xj passing through 
X must have length at least 21(x) — di — dj. Hence dq{xi,Xj) is at least 2/=,, — di — dj, where l^ 
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is the maximal label on a geodesic path from Xi to xj. Thus dq{xi,Xj) is strictly larger than 
\di — dj\ if l^: > diV dj. Obviously this is the case since the Xj's are pairwise distinct and have 
neighboring vertices with strictly larger labels, by construction. 

Finally, we infer that dq(xj, xj) + di — dj G 2N from the fact that Xi can be colored according 
to the parity of di — dj. If Xi and Xj have same color, then they are at even distance and di — dj 
is even, and otherwise both d(^[xi,Xj) and di — dj are odd. □ 

Note that, due to property 4. of Lemma El the set of half-edges of m is identified with the 
set £"1/2 of edges of q with the special orientation of Sect. 12.21 in such a way that half-edges of 
m incident to /i(m) correspond to elements of £'j(q, x, D). This will be useful in the sequel. 

We now state the main result of this section. 

Theorem 4 Let k >\ and > 0. The mappings . and <!>° . are inverse of each other, and 
therefore are bijections between the sets Q° ^ and LM°^^. Moreover, if (m, [1]) = ^(q, x, D), 

L \F{q)\ = |E(q)|/2 = |E(m)|, more precisely \Ei{(i,x,D)\ = degin(/j(m)), 1 <i<k. 

2. |F(q)| = |y(m)|+fc, 

3. For every 1 < i < k and d> 1, 

\{e £ Ei{q,x,D) : (iq(e~, j;^) = d} 

(note that the last appearing quantity is independent of the choice 0/ 1 G [I]). 

Properties 1,2,3 are easy consequences of the construction and Lemma [5J They are far from 
being the only important properties of f^- Other consequences will be derived in Section [2.71 

Proof. It is easy to see why 'I'" j^i^g ^(m) [1])) = (rH) [1]) from the construction of the 4-valent 
faces of $° fc(m, [1]) depicted in the previous proof. We just notice that in each of the three cases 
that are detailed (e G E{m) is such that l(e") = l(e+), l(e") = l(e+) + 1 ox l(e~) = l(e+) - 1), 
the edge e that is deleted in the end is precisely re-added by the construction of ^ (the first 
case is that of a confluent face, and the two others correspond to simple faces). The labels 
remain unchanged in this step thanks to the last sentence of Lemma [2l 

We now prove that ^° ^(^g ^(q, x, D)) = (q, x, D). To see this we must check that the edges 
of q can all be obtained as edges between a corner and its successor in "if" ^(q, x, Z)) = (m, [1]). 
Necessarily, edges must join corners with labels differing by 1, and cannot cross the edges of m, 
so that they all lie inside some of its faces. Moreover, inside the face /j(m), all corners with 
label di + 1 must be re-linked to the vertex Xj. So let e, e' £ /j(m) be linked by an edge of q, 
for which 1(e) = l(e') -|- 1 = / -|- 1, but such that s{e) 7^ e' . We choose moreover e, e' in such a 
way that / is minimal, that is, all the edges between a corner with label < / and its successor 
are edges of q. In this case, the paths s{e), s^(e), . . . , Xj and e', s(e'), . . . ,Xi are paths in q, and 
cannot cross the edge between e and e'. On the other hand, s(e), s^(e), . . . must lie in [e,e'], 
while s(e'), s^(e'), . . . must all lie in [e',e]. In particular, adding the edge between e and e' to 
/j(m) breaks the latter in two connected regions, both of which have to contain Xi, and this is 
impossible. □ 



lee fi{m): 1(e) - min l(e') + 1 = d\ 

[ e'6/i(m) J 
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2.6 Rooting 

We let Qg^k be the set of 3-tuples (q,x, D), where q is an element of Qg, and if q° denotes the 
associated unrooted map, (q°,x, D) G Q^^- Similarly, we let LMg^^ be the set of pairs (m, [1]) 
where (m°, [1]) e LM°_fc. 

Let (q, X, D) E Q^.fc, and (m°, [1]) = ^.(q°, x, D). With root edge 6=,, of q, corresponds nat- 
urally one of the corners of (m, [1]), i.e. it is the unique corner such that the edge drawn between 
e and s(e) is either equal to e* or its reversal when performing the converse construction. We 
root m° at this corner, obtaining a rooted map m, and let ^'^^^(q, x, D) = (m, [1]). Notice that 
since any edge of q joins two vertices with labels that differ by 1, the two possible orientations of 
a base edge of q yield two distinct rooted and fc-pointed quadrangulations (q, x). The following 
corollary of Theorem U] follows. 

Corollary 1 The mapping ^ g^k '■ Qg,k — > LM^jt is two-to-one. 

Besides getting rid of symmetries, the nice feature of rooting is that it will allow us to choose 
a particular labeling function inside its class, as we will see in Sect. El 

2.7 Bounds on distances 

The construction of ^° ^ allows us to get useful bounds on distances between vertices in the 
quadr angulation, from the label function on the associated map. 

Lemma 3 Let (q,x,Z)) G Qg,fc ^^^^^ ("^i [1]) = ^g,fc(q) x, D). Take two half-edges ei,e2 G 
Ei{q,x, D) for some i. Still denoting by ei,e2 the half-edges o//i(m) associated with ei,e2 by 
^g^kj we have 

|l(ei) -1(62)1 < dq(e^,e^) < I(ei) + 1(62) -2 min 1(e) + 2 

eG[ei,e2] 

(here the interval [61,62] is the one in m). 

Note that we could have interchanged the roles of 61 and 62, so that we have in fact 

dq(ej", 6^) < 1(61) + 1(62) — 2 max ( min 1(e), min 1(6) I + 2 

\eG[ei,e2] ee[e2,ei] / 

Proof. The lower bound is obvious from the fact that for e G -E.j(q, x, D), 1(e) — l(xj) = dq{e~ ,Xi) 
and by the triangular inequality. For the upper bound, consider the last half-edge e, among the 
sequence s"^{ei),m > 0, that lies in [61,62]. Necessarily, 1(6) < 1(62) as otherwise, the successor 
of 6 would be in [6,62]. From this, it follows that s(e) = s^^^^)~'(^)"'""'^(e2). Therefore, there exists 
a chain in q from e^ to e^ with length (l(ei) — 1(6)) + 1 + (1(62) — 1(e) + 1), which yields the 
result. □ 

3 Decomposition of labeled maps 

Our next task is to make one more step in our bijective study, by pulling apart elementary 
pieces that compose an element of LM^^yfc. Our presentation resembles that of [37], the two 
differences being that our maps are rooted, which allows us to get rid of the symmetry factors, 
and, more importantly, the vertices are labeled. In the following, it will often be implicit that 
(5) ^ {(0) 1)) (0) 2)}, which need special treatment. 
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3.1 Map reductions 

We let Mg be the set of rooted maps of genus g with k indexed faces. Elements of ^ for 
fixed k but with a large number of vertices, show a dendritic pattern branching on a graph with 
a simple structure. In order to study this in a precise way, we are going to perform several 
reductions on the maps. We let M-^ and M-^ be the subsets of ^ consisting of the maps 

^2 

in which all vertices have degree at least 2, resp. at least 3. It is easy to see that M-, is an 

infinite set (take a simple loop with an arbitrary number of edges). On the other hand, is 
finite according to the combination of Euler's formula 

\V\-\'E\ + k = 2-2g (6) 

and the fact that 

2|E| = ^deg^;>3|y|, (7) 

which entails |E| < 3k + 6g — 6. Note that equality holds if and only if every vertex has degree 3 
exactly. We let ^ be the subset of M-^ for which this extra constraint is satisfied. According 

to these definitions, the set should be empty. We make an exception by letting Mq 2 = 
consist in the two loop maps made of a rooted loop with an origin, bounding two faces with 
degree 1 (there are two such maps according to whether the face indexed 1 stands to the left or 
to the right of the root edge). Elements of will usually be denoted by the letter m, and 
their root by e,,,. 

There are natural projection maps 

which roughly consist in taking the largest subgraph belonging to (resp. M^^), for some 

appropriate definition of subgraphs, where merging edges are allowed. More precisely, starting 
from (a graphical realization of) m G we delete one by one the edges incident to vertices 

with degree 1 until no such vertex remains. This operation does not disconnect the graph, nor 
does it modify its genus or the number of its faces. The outcome is thus a map m-^ with k 
faces and without vertices of degree 1, and which does not depend on the order of the edges that 
were removecH. If m has root e^,, we naturally root m-^ at the first element of ipl^{e),i > 
that belongs to i?(m-^). This yields an element m-^ G ■'^g fc" 

Likewise, starting from m G ^g,ky we can remove all the vertices of degree 2 in m-^: we 
concatenate the half-edges ei, . . . , of a simple chain linking two vertices of m-^ of degrees 
> 3, and traversing only degree-2 vertices, into a single edge e = ei . . . e,.. We root the resulting 
graph at the one edge = ei . . . Sr in which the root of m-^ is involved. This does not modify 
the genus nor the number of faces of the map, and these merging operations commute, and thus 
yield a well-defined element m-^ of M^^. It is also obvious by definition that (m-^)-'^ = m-'^. 
There are natural inclusions V{m-"^) C y(m-^) C V{ui) and £'(m-^) C E{m). The map m-'^ 
is sometimes called the homotopy type of m. See Figure [5l 

■^Alternatively, elements of -E(m-^) are those that can start an infinite non-backtracking chain ei, 62, • . ., i.e. 
such that ei 7^ e^+i for i > 1. 
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Figure 5: Top-left: a portion of a map m, emphasizing the tree components branching on 
m-^ depicted on top-right. Bottom: the projection m-^. The picture shows how the various 
projections are re-rooted depending on the location of the root of m. 



3.2 Some notations for continuous paths 

Certain continuous functions will be crucial to encode maps, so let us introduce some notations. 
A path is a continuous function {w{s),0 < s < t{'w)) taking its values in R, such that w{0) = 0. 
The number t{w) > is called its duration, and w = w{t{w)) its terminal value. We also let 
w = w{t{w) — ■) — w he the reversed path (beware that this does not mean the same as the 
reversal of an edge in a graph). We let W be the set of paths, and make it a Polish space by 
endowing it with the distance 

dist(i(;, w') = I sup \w{s A t{w)) — w'{s A t{w'))\ I V \t{w) — t{w')\ . 

\s>0 J 



w^s) if < s < r(w) 

w + w'^s — t{w)) if t{w) < s < t{w) + t{'w') = t{ww') 



If w, w' G yV, we let ww' be the concatenation of w and w' defined by 
ww'{s) = 

For {g, k) ^ {(0, 1), (0, 2)} introduce the sets 

Cg,k = {(m, M*) : m G Mf„ [ = (f, e G E{m)) G < t* < r^"] 

endowed with the natural product topology. We make the exceptions Co,i = {(c, z) G W^} and 
Co,2 = {(m, l,t*) :m G M^^, [ = {{c, z),l'' ,(^*) G x x >V^ < < r(c)} . 
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For (^f, k) ^ {(0, 1), (0, 2)} the generic notations are 

r = K,c^z^), re = rK), r^ = r(c^), n= J] T^ (8) 

For g = and /c = 1, 2 we let r = t(c), for k = 2 we also let rj = r + r^* + r^* . Our next task is 
to bijectively encode elements m of LM^^^ with certain elements of Cg^k, whose first component 
is the homotopy type of m. 

3.3 Level sets of Mg^k ^ M|| 

Let us see how to get back from the projection m = m-^ to the initial map. First of all, every 
e G E{m) must be split into a certain number of edges ei . . . e^^ , with r,. > 1, to yield the element 
m-^, and obviously = rj := re (where e = {e,e}) so what is needed to do this are positive 
integers (re,e G E(m)). 

Then, consider elements e G E{m) such that the first element of — belonging to 

i?(m-^) is one of the edges ei, . . . ,6^^. These edges form an interval, and we let e(0) be the 
first element in this interval, and label the others as e(i) = ^mit{0)), < i < — 1. With this 
labeling, we see that every half-edge of m is of the form e{i) for a unique e G m. Moreover, 
the half edges of the intervals (ej_i,ej),l < j < Ve must form an acyclic graph, i.e. a tree, 
possibly reduced to a single root vertex (in particular, their set has even cardinality and is 
stable under edge-reversal). Here, we take the convention that (eo,ei) are the elements of the 
list e(0), e(l), . . . that precede ei strictly. Otherwise said, we can interpret e(0), . . . , e(r^ — 1) as 
a forest of trees, with roots bound together by the floor half-edges ei, . . . , e^j,. If is the 
root edge of m-^, then it deserves special attention. Indeed, in this case, one of the edges of the 
interval [e*(0), e* (r^* — 1)] must be distinguished as the initial root edge of m, which amounts 
to distinguishing an integer t* G [0, r"^). 

The contour function associated with this forest is defined by 

c\i) = dM^~.V{m^^)) - c\i) , < i < T% 

where 

(fii) = \{j: e, G [e(0), e(i)]}| = - inf .c^(j) , < i < - 1 , 

0<j<« 

and c^(r^) = r^, is the number of trees in the forest already explored before time i. We extend 
into continuous functions (still denoted by c^,c^) by linear interpolation between values 
taken at integers. The result is a concatenation of Dyck path^ encoding the trees of the 
forest, and achieving a new negative infimum when a floor edge of the forest is traversed (note 
that the terminal value —r^ of is attained for the first time at the final time r^). It is an easy 
and well-known fact that the contour function entirely determines the forest. In this way we 
have shown 

Lemma 4 Assume {g,k) ^ {(0, 1), (0, 2)}. An element m G M^^fc is unambiguously described 
by its projection m = m-'^, a family (c^, e G E{m)) G yV^^'^') such that is the concatenation of 
Te Dyck paths, and an integer t* G [0, t^* ) . 



*Here we will call Dyck path a piecewise linear element c £ W starting at 0, with slope ±1 between integer 
times, and which attains the value —1 for the first time at r(c). 
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Figure 6: The contour process c^* associated with the edge e* in the reconstruction of the map 
m from m-'^ as depicted in Figure O The process —c^* is drawn in thick dashed hues. The 
location of the root of m is distinguished by the time t*. 



See Figured] for an illustration. The cases {g, k) = (0, 1) or (0, 2) are special and requires a 
little discussion. An element of Mo,i is simply a rooted planar tree with at least one edge, and is 
described by a single Dyck path c with duration r(c) > 1. When (^f, k) = (0, 2), the construction 
of m-^ from an element m G Mo,2 shows that the root vertex of m-^ can be any of the vertices 
of the loop constituting m-^. This difficulty does not appear for other values of {g,k), since 
the root edge of m-^ has an unambiguous origin in the original map m. While this does not 
affect the final outcome considered as a map, there are several non-equivalent ways to obtain 
the same result. To lift the ambiguity, we use the extra convention that the root vertex in m-^ 
is e~ , where e is the first element of S(m-^) in the sequence (</'m(6*))* ^ 0). Otherwise said, 
we single out the tree component that contains the root and branching on the unique loop of m, 
by saying that its root vertex must be the root vertex of m-^. The other tree components are 
obtained by taking concatenations of Dyck paths so that an element m S Mo,2 is unambiguously 
described by the element m-^, a Dyck path c, an integer time in [0, t(c)) and two paths c^* , cf", 
the first being a concatenation of r^.^ — 1 Dyck paths, and the second the concatenation of r^, 
Dyck paths. 

3.4 Labeling the maps 

In order to complete the description of elements of LM^^/t, we need to incorporate the labeling 
function [1]. We can use the origin of the root edge e* of m = m-'^ as a reference for [1], i.e. 
we can choose the element 1 G [1] such that l(e~) = 0. This convention will be in force from 
this section onwards. It is then natural to introduce a labeling function T for every e G E(m), 
defined by 

lt = l{t{i))-im), 0<i<r% (9) 

and l^c = l(e(T^ — 1)"*") — l(e(0)), where the notation e(0), e(l), ... is as above the list of edges 
of the trees branching on e. We interpolate linearly between integer times, in order to define 
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a continuous piecewise linear function with slope in {—1,0,1}, starting at and with same 
duration as c^. Another process that will be of crucial importance is the process 

k = K^inMO))) < i < deg^(A) - 1 , (10) 

where /* is the face of m incident to e*, and l^^g (j^) = 0. As usual, it is turned into an element 
of W by linear interpolation. In the case where m has only one face, we see that / provides 
a natural exploration of the labels of all vertices of m in the fatgraph face order. It is easily 
checked that I is the concatenation 

i = i^^i^^M__j^i''-'^*'-'M^ (11) 

where f* is the face of m incident to e^=. In order to study scaling limits, we will still need to 
decompose T into two more elementary pieces, which is the purpose of the two next sections. 

3.4.1 Walk networks 

A walk is an element w G W taking integer values on interpolating linearly between these 
values, and satisfying \w{i) — w{i — 1)| < 1 for 1 < i < r := t{w) (where it is assumed that 
r > 1). Let W be the set of walks. A walk from to x G Z is a walk with final value w = 
their set is written Wq^o;. 

Let m G for some g,k > 0. We let £{V) and £{E) be the spaces of functions ^(m) R 
and E{m) M respectively. We also let ^o(^) be those functions 1 of £{V) with l(e~) = 0. A 
function G £{E) is antisymmetric if 6{z) = —0(e) for every e, and we let i-{E) be the space of 
antisymmetric functions. For 1 G ^iV), we let VI G £-{E) be defined by Vl(c) = l(e+) — l(e~). 
As m is connected, the kernel of the linear map V is the space of constant functions, so the 
restriction of V to io(^) is an isomorphism of vector spaces. A function 9 in the image of 
V is called a potential. 

A walk network on m is a collection {w^, e G E{m)) G W^^™) such that 

w'it) = w'{re -t)-w', < t < re = t{w') = r(w^) , 

and such that the antisymmetric function 9 = w = (w^, e G E{m)) is a potentiajl. We let WNm 
be the set of walk networks on m. If we are given an antisymmetric function G {E) assuming 
integers values, we let WNm(0) be the set of walk networks (u;^, e G E[m)) G WNm such that 
w = 9. 

If (m, [1]) G LMg^fc, and using the notations of Section [THl we set, for every e G Eixn) of the 
form e = ei . . . Cr-g, 

w\i) = \{eT^^)-\{e^), 0<i<re-l, 

and w^{re) = l(e^) — l(e^). This defines a walk with duration r^- Moreover, the family 
{w'^, e G E{m)) is a walk network. The reversal property tL'^(-) = w^{re — ■) — {v^ is clear since 
e = Cr-g . . . ei, and by definition the function {iv^, z G E{m)) is VI|v{m)- 

'''Alternatively, walk networks on m are in correspondence with labeled maps m £ such that m-"^ = tn. 
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3.4.2 Discrete snakes 

Finally, we say that the pair (c, z) £ is a discrete snake on r > 1 tree components if 

• c is a concatenation of r Dyck paths, 

• z is taking integer values on Z+, interpolating linearly between these values, with inte- 
ger steps in { — 1,0,1}, satisfies t{z) = t(c), and z{i) = z{j) whenever c(i) = c(j) = 
minjAj<fc<ivj c(/c) (which means that i and j are times where the same vertex of the forest 
encoded by c is explored), 

• moreover, z{i) = whenever i is an integer time at which c achieves a new minimum (this 
means that i is a time at which a new tree component is explored, and the label of the 
root is fixed to be 0). 

In words, a discrete snake is the process of heights and labels when exploring an integer-labeled 
forest in contour order, where the labels differ by at most 1 in absolute value between adjacent 
vertices, and the labels of the forest's roots is 0. We let Sr be the set of discrete snakes on r 
tree components. 

If (m, [1]) G LMg we set, whenever < i < rMs such that e(i) E (em_i,em] (recall the 
notations of Sect. IS.Sp . 

z^i) = I(e(i)) - l(e^) , 

and we interpolate linearly between integer times. The pair {c^, z^) is then the discrete snake 
associated with the forest branching on the edge e, where the labels in each tree component have 
been shifted so that the roots have labels 0. With the data {w^,c^,z^) defined so far, the label 
function introduced at the beginning of the subsection is recovered as 

l\i) = z'{i)+w\(f{i)). (12) 

We also use to define the function T out of any (m, I, t*) G C^^/t, where (f{t) = — info<s<t c^{s) 
for < s < r^, and define the function / as in (jlip . 
Summing up the study of this section, we obtain 

Proposition 1 (i) Let {g,k) ^ {(0, 1), (0, 2)}. Then LMg^^ is in bijection with the subset 
C^^^ C Cg^k of elements (m, l,t*) such that 

1. {w\ee E{m)) G WN^, 

2. (c^2;^) G Sr^ for every e G E{xn), and 

3. t* is n integer in [0,r(c^*)). 

(ii) LMo,i is in bijection with the set C^^^ = {(c, z) G Si : t(c) > 1}. 

(iii) LMo,2 is in bijection with the subset C^2^ C Co, 2 of elements (m, (, t*) with 

1. w^* G Wo_>o and w^* =w^*, 

2. (c, z) G Si and t* G [0, r(c)) is an integer, 

3. {c'*,z'*) G Sr,,^i and {c'%z'*) G S,.^^ . 
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4 Discrete and continuous Boltzmann measures 

We now explore natural measures and probability distributions on tessellated genus-*/ maps, and 
their continuous counterparts. 

We define the measure Q^^fc on the set Qg^k as the cr-finite measure assigning mass 12~l^('')l 
to every element (q, x, D). Otherwise said, if we recall the notations of Section [1.41 one has 

Q,,fc(d(q, X, D)) = Q,(dq) (yqMq)®'=(dx) #v(n,^){dD) , (13) 

where #x)(q,x)(d-C') is the counting measure on the (possibly empty) set of delays in q between 
sources x. Up to Sect. 14.41 we assume {g, k) ^ {(0, 1), (0, 2)}. 

According to Theorem [J] and Corollary [H the push-forward of Q^^^ by the mapping ^g^k is 
the measure on the set LMg^^ that assigns mass 2 • 12~l-^^'"^l to any element (m, [1]). Using the 
considerations of Section [3l we see that this measure can in turn be obtained out of measures 
on the more elementary pieces that compose a labeled maps (walk networks and snakes). More 
precisely, consider the bijective mapping '^g^k '■ LM^.fc — > C™^^ that is described in Proposition 
III Define 

:= Sg,fc o ^iig f, : C^l^ ■ 

Then the image measure LMg^^ of Qg^k by ^g^k assigns weight 2 • 12""^'/^ to (m, (, t*) G C™^^, 
as one checks that T(/2 = |E(m)| = |-F(q)| = V"q — xid) whenever (q,x, D), (m, [1]), (m, (, t*) are 
associated by our bijections. We rewrite this measure in the form 

LM,,fc = 2 LMn, (14) 

where 

LM„ = 5|„|WN^(d(^^),es(„)) (g) DS,(^c)(d(c^z^))#,(,c,)(d^*) 

and 

• WNm is the measure assigning weight 3~^eGE(m)'"e ^ walk network {w^) £ WN^^, 

• DSr is the measure assigning weight 2~'^'''^-*3~^'^'''^-*~^-*/^ to an element {c,z) G S,., 

• lastly, #t{c'*) is the counting measure on {0, . . . ,r(c^*) — 1}. 

Individually, these measures are c-finite measures, which are finite except for WNn^. The 
measure DS,. is a probability measure, under which the r trees involved in the forest are inde- 
pendent Galton- Watson trees with geometric (1/2) offspring distribution, and with uniformly 
distributed labels among the allowed ones, conditionally on the forest. To see this, recall that 
(c, z) G Sr encodes a forest with r tree components, and t(c) oriented edges (counting the r floor 
edges between tree components). A component of this forest, which is a labeled tree with m' 
edges say, contributes a probability factor 2~^"^ -ig-m weight 2~'^*^'^)3~^'^^'^'*~'')/^. Also, 

the measure DSr(d(c, /;))#^(c) (dt*) is an infinite measure. 

4.1 Walk networks 

For simplicity we let in this section V = V{m),E = £'(m),E = E(Tn) for some m G M^^. We 
let WNni be the image measure of WNn^ under {w^, t £ E) ((''e, e G E), {w^, t G E)). 
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4.1.1 Gaussian laws on potentials 

Fix m G fc, and let = t~ be the root vertex. Recall the notations of Section 13.41 We 
endow ^o(^) with the scalar product (x, x') = Ylv^v* x^x'^, and endow the vector space Vy := 
V^o(^) C (--{E) of potentials on m with the scalar product 

(Vx, Vx') = (x, x') . 

This space is isometric to ^o(^)) and of dimension — 1 = dim£o(V^)- We let Ay be the 
measure on Vy which is the image of riijeFVfi;*} ^ • ^o(^) ^-{^)^ a-^id call it the 

Lebesgue measure on Vy: it assigns mass 1 to the unit cube delimited by the vectors of the 
orthonormal basis (Vlt,,f / v*). For (re,e G E) a family of positive weights, we also let 

'S'r(x) = ^ xl/r^ , 

eGE 

which defines a positive definite quadratic form on Vy. The Gaussian distribution on Vy 
associated with Sj. is the measure 

AAr-^ldx) = ]— exp (-^^) yd^Av(dx) , 

(2^) 2 V 2 ; 

which is a probability distribution. One should be careful that det £r is the determinant of the 
quadratic form defined on (Vy, (•,•)) and not the whole space i-{E). 

Lemma 5 Let ST(m) denote the set of spanning trees ofm. Then, one has 

det£,= n (^-)"'- 

aeST(m) eGE(a) 

Otherwise said, this determinant is a partition function for spanning trees on m with weight 1 jr^ 
for the edge e. 

Proof. In the orthonormal basis (Vl{^},u ^ v^) of (Vy, (•,•)), it is easy to check that the 
matrix {niuv, u,v £ V\{v*}) of the quadratic form Sj. is given by m^^ = X]e-e-=v ^/^^ ^^^^ ~ 
~X^e-e-=u e+=v l/'^e- Otherwise said, this is a minor (excluding row and column corresponding 
to u*) of the weighted Laplacian matrix on V with weight l/r^ on edge e, and one concludes 
by the matrix-tree theorem (and the result does not depend on the root of m but only on the 
map). □ 

The following result will also be useful when estimating integrals with respect to Ay- For a 
spanning tree a S ST(m), we privilege the orientation of its edges Ei/2i.'^) by letting them point 
away from the root vertex v^. If y £ R^i/^^"), we define a function 6y G ^~{E) extending y by 
letting ^y(e) = —y^ for e G £'1/2(1) and then 

where the sum is over the half-edges e' G E[a) of the unique oriented simple chain going from 
e~ to e"*". It is easy to see that 6y G Vy is a potential, and any potential is recovered from its 
restriction to £1/2(1) by applying 9°-. 
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Lemma 6 Let H : R^ — > be measurable. Then for a G ST(m), 

/ Av(dx)/7(x)= / ^dyH{9-). 

Proof. By definition 

[ Av(dx)i/(x) = [ dzH{Vz) . 

The result follows from the change of variables y(e) = Vz(e), e G £^1/2(1)1 provided we can justify 
that its Jacobian is 1. This results from the fact that this change of variables is linear, and that 
its representative matrix in an appropriate basis is triangular with only I's on the diagonal. 
To see this, perform the search-depth exploration of a starting from any half-edge ei G E{a) 
with ej~ = f*, and letting, for i > 2, = 99^^*^ (ei_i) where m(i) is the first integer j > 
such that (93-^(ej_i))"*" is a vertex that does not belong to the set {v^, , . . . , of already 
explored vertices: this is possible as long as i < \V\ — 1. This gives a labeling ei, . . . , en/|-i of 
the half-edges of a pointing out from u^,, such that the source e~ is always a vertex belonging to 
the set {e^ : j < i}, while the target is a newly explored vertex. □ 



4.1.2 Key lemmas 

Define the quantity 



E n 2vrre, rG 
oeST{m) e^E{a) 



which is a partition function for the complement of a spanning tree, with weight r^ on edge 
e, and defines a homogeneous function of degree |E| — \V\ + 1. We define a measure PN^^ on 
X Vy by 

PN™(dr,dx) = ^i=AA-(dx). 
Let Ae = {u = (tie, e G E) G : X] 

egE ■"e — 1} be the |E| — 1-dimensional simplex indexed 
by E, and A a denote the (scaled) Lebesgue measure on Ae, satisfying 



/ H{r)dr= [ pl^l-Mp / AA(du)if(pu) 

7r5 Jr. Ja-e 



for every integrable H. In particular, taking H{r) = exp(— ^^^e '''e), we see that its total mass 
equals 1/(|E| - 1)!. 

We define the measure RnJ^^"* on /oAe x Vy as the image measure by (u, x) 1-^ (pu,x) of 



pl^l-^^^^^AC(dx), 
hence giving a disintegration of PN^^ with respect to p{r) := X^eeE'^e 
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Lemma 7 For u G ^Ae, we let u be the point o/nAEPlZ^ that is closest to u for the Euclidean 
metric (taking some arbitrary convention for ties). Then for any bounded function H : M, 



< (sup|i?|)e„ 



where the sequence {en,ri > 1) is independent of H and goes to as n ^ oo. 
Proof. To see this, note that 



XA{du)H{f^)= J2 ^(v)Aa({uG Ae : ?^Ti = v}) 

vGnAEnZE 



(15) 



Now, as long as v G kAe H has all its coordinates positive, i.e. does not belong to the 
boundary ndA^, and therefore stands at distance at least 1 from it as it has integer coordinates, 
the sets {u G Ae : nu = v} do not intersect c?Ae, and have equal AA-nieasure by an obvious 
property of translation-invariance. There are (|e|1^i) elements in nAEHN^, which is the number 
of integer compositions of n into |E| positive parts. We conclude that for any v G tt-Ae H N^, 



< Aa(Ae) 



n- 1 
El - 1 



Aa({u G Ae : nu = v}) < Aa(^i/„(5Ae)) , 



where Vi/„(9Ae) is the l/?i-neighborhood of (9Ae, and AA(^i/n(<9^E)) ^ as n ^ oo. Since 

(|E|"_\)~-''^'"V(|E|-l)=nlE|-^AA(AE), 

we conclude that Aa({u G Ae : nu = v}) = (1 + o{l))/n''^^~^ as n ^ oo independently of 
V G tiAe n N^, and since \ndA-E, D Z^| = o(nl^l~^), the conclusion follows easily from ()15p . □ 

Lemma 8 Set ? = (8/9)^/^. Let Ga,a > be a uniformly bounded family of functions on 
X such that 

Ga(ra,Xa) ^ ^(r, x) 

whenever Va/ {2aY/'^ r and Xa/(?a^/^) ^ x G M.^ . Let also pa/(2a)^/^ p with pa G N. Then 

k+2g-l _ ( ) 

2 2 O 4 

Proof. Let Wg^^. be the set of walks with duration r that end at x, and Pn{N) = S^'^'IWq^^^I 
be the probability that a uniform walk of duration r ends at iV. Note that A;+2y— 1 = |E| — |y| + l 
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by ([6]), and write the left-hand side in the statement as 

|E|-|V^|+l 

|E|-|y|+l 

E E G«(r,Vy) n Pn(Vy(e)) 



|E|-1 |\^i + |E|-3 
2 4 rGpaAEnNEygZ^Xf"*} eG^;!/. 



~ Pa 
a— >oo 



a— >oo 



AA(du) / dyGj/^,VLaVVj) X{<^a'"ppr^i.){y[a"\y\{t)), 



where we made use of Lemma [7] for n = pa at the third step of the computation. Here, [xj is 
the vector whose components are the integer parts of the components of x, and we write u(e) 
for the coordinate of u corresponding to the edge e = {e, e}. 

By the local limit theorem |39l Theorem VII. 3. 16], there exist^ some C G (0, oo) such that 
for every n > 1, G Z and letting x = N/{n'^/^J2/3), 



1 + X 



2 ■ 



Prom this, and the easily checked fact that ndA^ is at distance Ce > (independent on n) of 
the set {nu : u S Ae : nu ^ ndA^} (so that the components nu^ are bounded away from for 
nu in this set), it follows that 



?«'/'pp';n(e)(VLa'/\yJ(e)) < Cmin 



where = |Vy(e)|. Combined with the fact that the Ga are uniformly bounded allows us to 
bound the last displayed double integral by 



Cmax(p3|E|/2-i^^|E|/2-i) /■ Av(dx) / AA(du)n f ^^KM} + -^^K 
Developing the product and integrating in u, we finally obtain the bound 

Cmax(p3|E|/2-l^p|E|/2-l^ f -Q ^in(^-2^i). 

Taking a spanning tree a € ST(m) and using Lemma O this is bounded by 

max(pW-i^^|E|/2-i) / ^^^^^j" ^ ,nin(x7^1)) , 



<1} 



(16) 



^Until the end of the proof, the same notation C will be used for different such bounding constants. 
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and this is finite. Therefore, we can apply dominated convergence in the integral expression 
above and pass to the limit as a ^ oo to get 

|E|-|V|+1 

WNn,(G,l,(r)=pJ 



|E|-1 |V-| + |E|-3 

2 2 a 4 



J^l-^ / AA(du) / Av(dx)G(pu,x) n 



cG-Bi/2 



-x^/{2pu,) 



exp ' — ^ ^ ' 



r r exp i — i 

p|E|-W AA(du) / Av(dx)— A_^=^G(pu,x), 



and the result follows from the definition of A/"^ and Z™{r) and from Lemma O □ 
Lemma 9 Let Ga,a > be functions x — > M such that 



|Ga(r,x)| <Kexp 



re 

eeE 



for some K,c > 0, and satisfying the same convergence assumptions as in the previous lemma. 
Then 

2 2a 4 a^oo 

Proof. This is an immediate consequence of the previous lemma, using easy dominated con- 
vergence arguments along with the p-dependent bound (fT6l) . Details are left to the reader. 

□ 



4.1.3 Path network measure 

We are now able to describe the scaling limit of the measure WNm- Let PQ^^{dw) be the law 
of the standard Brownian bridge with duration t from to x, and PQ_^^{dwdw') be the image 
measure of PQ_^^{dw) under w ^ {w,w) S W^. 

Take an orientation £'1/2 of E. We define a measure by 

PN^{d{w',ceE))= /"pN4dr,dx) J] P^l^^^{dwW) , 

which is a measure on paths networks, i.e. collections {w^, t (z E) of paths with durations 
Te := t{w^) such that = and (we, e £ E) £ Vy is a potential. This measure does not 
depend on the choice of E1/2 because of the well-known invariance property [H] of Pq^^ under 
w ^w. 

We define a scaling operator (j)^ on W by 

/ q \ 1/4 

(j)<^(w) = ij-j It; ((2a) , 0<s<T{w)/i2a)^/\ 

We extend the scaling function 0^* to a function on powers of W in the natural way, i.e. it acts 
on every component by applying the scaling cj)"'. 
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Proposition 2 Let H be a continuous function — > M such that 

\H{{w',ee E))\ < Kexp j -c^r(u;') 

for some K,c > 0. Then 

k+2g-l 

IE ,V| + |E|-l WN„,(go</,") PN^(if). 

2— a 4 

Proof. We rely on Lemma [9] and the following statement. 



Lemma 10 For n G N and x G Z, let Po_^j. be the uniform distribution on Wg^^.. Let r > 
0, X G M 6e fixed and assume that ra, Xa are integers such that raj (2a)^l'^ — > r and (^l'^>d)^l'^Xa — >■ 
X. Then the 'probability measures 4''*^^^-^xa converge in distribution to PgL^j, in W. 

Proof. This is a classical elaboration on Donsker's invariance principle. See for instance [1'6\ 
Lemmas 4.5-6] for the scheme of the proof (in this reference, the author considers more general 
stable processes than Brownian motion but is also concerned only in the case where x = 0, 
however the proof in the general x case is the same). □ 

To end the proof of Proposition [2l we apply Lemma [9] to the function 
G„(r,x)= / (g) Pl'^.^^{dwW)Hicj)'^{w\tGE)), 

where as above the measure PQ^^{dwdw') is the law of {w,w) under Pq^^.. □ 

Remark. The path network measure PN^ actually has the following interpretation, which 
we state without proof as this is not needed in the present article. Assume that the numbers 
Te > 0, e G E are given and consider a graph with same structure as m but where the edge e 
is isometric to a real segment of length r^- This defines a metric graph [9l Chapter 3.2.2] rrir- 
Let Onxj, be the set of continuous functions on rrir that vanish at f*. With any path network 
(w^, t £ E) we may associate an element of in the following way: for x G rrir on the half-edge 
e and at distance < s < r^ from its origin, pick a chain ei, . . . , of half edges starting at 
and ending at e~ . Then let H(x) = ^ w^^ + w^{s). The reader will convince himself that 
this definition does not depend on the choice of the orientation of e nor on that of the chain 
ei, . . . , Cm, because of the properties of path networks. Thus, the measure 

A/7(dx) n ^'o^x.(dw^'d«;^) 

induces a law GFFm^(d-ff) on @mr- This law is what is referred to as the Gaussian free field on 
mr- We refer to |33] for a nice introduction to this object (here the metric graphs are really 
1-dimensional objects, so the associated free fields are well-defined random functions that are 
much less elaborate than higher-dimensional free fields on which [13] focuses). To explain what 
GFF is, imagine that we run a "Brownian motion" on rrir. It behaves as a linear Brownian motion 
when moving on an edge of rrir, and when it encounters a vertex of the graph, its excursions 
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out of the vertex choose independently at random in which incident edge they will diffuse. The 
diffusion is killed when it first encounters the vertex We let G{x, y) be the densit}{l| at vertex 
y of the occupation measure of the diffusion starting at x. Then GFFt^^ is the law of the centered 
Gaussian field {H{x),x e rrir) with covariance 

Cov {H{x), H{y)) = G{x,y) , x,yemr. 
4.2 Discrete and continuous snakes 

We now discuss the scaling limits of the measures DS^ on discrete snakes. This is more classical 
from a probabilistic point of view as we are back into the realm of probability measures. We 
define the continuous snake measures S,., r > 0, as the law of the pair ((^s)o<<j</-ii ('S's)o<s</"Oi 
where 

• [Bs, s > 0) is a standard Brownian motion in R, 

• {Is = inf o<s'<s i?s') s > 0), is the infimum process of B and I^^ is its right-continuous 
inverse, 

• conditionally on B, the process S is centered Gaussian in R, with covariance 

Cov {Ss, Ss') = inf {Bs" - Is") ■ 

sAs'<s"<s\/s' 

As an application of Kolmogorov's continuity criterion and the standard fact that trajectories 
of a Brownian motion are a.s. locally Holder continuous, a continuous version of S exists, and 
St is a distribution on W^. 

We rely on a version of invar iance principles such as [111 \3T\ [23], but holding for forests. It 
is easily derived from Theorem 3 of [35] by using the methods of [14[ Section 2.4] to derive the 
results for the contour process of the tree rather than the height process. There is a difference 
between the label processes used in this reference and the present paper (due to a different way 
of exploring the trees), but the proof are adapted in a straightforward way. 

We define the scaling functions , i/^" on W by 

^"(c) = j^^c{2as) , 0<s< T{c)/2a , 

and 



il;'='{z) = z{2as), 0<s<T{z)/2a. 

Lemma 11 Fix t > and assuTJic Td G M ^5 such that Vdj {2(i)^^'^ — > v. Then, the tHjEcl- 
sure (t?"^, 'i/'")*DSra converges weakly to Sr as a ^ oo. If Un is the uniform distribution over 
{0, 1, . . . ,n — 1}, we also have that the measure defined by 

DS,„(d(c, z))Ur^,){dt)H{r{c), r{z),t/2a) , 

for positive measurable H , converges weakly towards the measure defined by 

dtli 



[sAd{c,z))[ ^-^^^^ff^H{c,z,t) 
J Jr+ ■t{c) 



^The function y ^-» G{x,y) can also be defined as the solution of a Poisson equation A„,^H = vanishing at 
«*, for the appropriate definition of the Laplacian Am^ on regular functions of Qm^- 
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4.3 Scaling limits 

In this section, we assume {g, k) ^ {(0, 1), (0, 2)}. We define a measure on Cg^k by 

CLM^,A: = 2 Yl CLM„ (17) 

and CLMm is the appropriate image measure obtained by rearranging terms of 
<5n,PN^(d(«;^),e£M) (g) S,, (d(c^ z^)) dt* l|o<t.<...} • 

eGS(m) 

Note that the exponent > 3 of the first term defining LM^ ^ is replaced by a 3 here. We define 
the scaling operator by 

dm, l,f) = {m,{{r{w'),r{c'),r{z')),te E{m))X/2a) . 

Theorem 5 Assume {g, k) ^ {(0, 1), (0,2)}. Let H be a continuous function on Cg such that 

\H\ < irr-^exp(-/3T[) 

for some K, P > 0. Then it holds that 

k+2g-~l 

^^-^IMg^kiH o r) CLMg,k{H) ■ 



233-2 



3fc Sg 
' 2 a 2 • 4 



Proof. First of all, we check that for every r G N^, 

' " (18) 



DSr, (^exp (-^^i)) < K'exp ^-c^ja"^/^ ^ re j . 



Indeed, under DS,.(d(c, z)), the terminal time r of c has same law as the first hitting time of — r 
by a simple random walk on the integers. Thus r[ has same law under (^gg^ DS^, as the first 
hitting time of —2 X^^g^ by the random walk. Standard computations show that 



so that 



1 /l 2 

DSi(s"('^))= 0<s<l, (19) 

s 

DSr, fexp (-^Tij) = exp | -2(/> (-j X] | , 

where (f>{P) = ~ log(e^(l — \/l — e"^^)) is the Laplace exponent of the first hitting time of —1 
by the walk, so (f){f3/a) ~ (2/3/a)"'^/^ as a — > oo, hence (fT8|) . We can thus write LMyniH) as 

I WN4d(u;^)) j (g)DS,,(d(cSz^)) j Ur^.{dt*)r'*H{m, l,t*) , 



eS-E 
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with the notation C/„ of Lemma [TTJ Set 

G,(r,x)= (g) P5L.,,(g)DS, 



Because of the domination hypothesis on H, the fact r^* < T(, and inequahty ()18|) . the function 
satisfies the domination hypothesis of Lemma O Applying the latter and Lemmas (TU] and [TT] 
allows us to conclude that for every m G fc' 

(2a)-^ iK|' !E|+iv - '-Mm(i^°g") - CLM„(i7) 
2 2a 4 

Note that for m ranging in M|^, the quantity |1^| + |E| appearing in the scaling factor is maximal 
if and only if m belongs to the set ^ of trivalent maps of genus g with k faces. Indeed, 
|E| — \V\ = A; + 2^ — 2 is a fixed quantity, so that \V\ + |E| will be maximal if and only if |E| is 
maximal subject to the constraints. On the other hand, we already noticed that |E| < 3k+6g—6, 
with equality if and only if m G M^^^, in which case {\V\ + |E| +3)/4 = 5g/2 + {5k - 7)/4. Since 

M^^ is a finite set, we deduce that the only terms in the sum (jl4p remaining in the scaling limit 
will be those belonging to ^. Hence the result. □ 

Remark: the Hurwitz measure. It should be noted that a non-labeled analog of CLMg^^ has 
already been considered in the literature \37\ I38|. in problems dealing with the enumeration of 
maps, random matrices, and related algebraic geometry problems. Specifically, the first marginal 
Pf of the measure St can be considered as the law of the contour process of a random forest. 
The measure considered in the above references is 

where the first sum is over unrooted maps in ^. (but still with indexed faces) and the factor 
involving Aut (m), the automorphisms of m, allows one to get rid of symmetries of the graph. 
Note that we can regroup the factors in the above measure as (^^g^^^^ dr^Pr^{dc^)Pr^{dc^). In 
this form, and by using an interpretation in terms of trees of a well-known decomposition of the 
Brownian excursion due to Bismut (see Le Gall [25]), we see that we have replaced each edge m 
by a doubly-marked Brownian continuum random tree taken under the Ito measure, using the 
two marks as endpoints of the edge. The measure CM° ^ appears as a scaling limit of a random 
element of ^ where k remains fixed but the number of edges tends to infinity. Of particular 
importance is the Hurwitz measure, which is the image of CM° by the map 

(c^ c G S) 1 — ^ I ^(^')' ^<i<k 

>ee/i(m) 



as its Laplace transform has an interpretation in terms of generating series related to intersection 
theory on the moduli space of curves [55] . 
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4.4 Exceptional cases 

This section is devoted to the study of the planar cases {g, k) G {(0, 1), (0, 2)}. An element (m, [1]) 
of LMo,i is just a rooted labeled tree and is described by a snake Ho,i(m, [1]) = (c, z) G Si. The 
image measure of Qo,i by ^o,i •= ^o,i ° ^o,i will assign mass 2 • 12""' = 4 • 2~2"~i3~" to each 
discrete snake {c,z) with duration r(c) = 2n + 1, coding a labeled tree with n edges. This is 
precisely the measure 

LMo,i :=4DSi. 

We define its scaling limit CLMo,i as follows. Let N{dc) be the Ito measure on W, governing the 
intensity of the excursions of standard reflected Brownian motion (with normalization N(t(c) E 
dt) = (27rt^)~^/2dt). Let N(d(c, z)) be the measure on such that under N, c has "law" A^, 
and conditionally on c, 2: is a centered Gaussian process with covariance 

Cov (z(s),z(s')) = inf c(s") . 

sAs'<s"<sVs' 

Finally, set CLMo,i = 4N. 

Proposition 3 Let H be a continuous functional on such that 

\H{c,z)\ < AVexp(-/3r), 
where t := t(c), for some K,f3 > 0. Then 

{2a)'/\MoAHo{'^^,^|J^)) CLMo,i(i/). 

Proof. First integrate out the label part, i.e. assume H is only a function of the first component 
and write 

{2a)'/\Mo,i{H{r{c))) = 2{2af'' f dt-l^Catp.^ Pj*l(i/(r (c))) 

where is the uniform distribution over Dyck paths with duration 2n+l, the latter set having 
cardinality Cat„, the n-th Catalan number. It is well-known (and can be obtained from Lemma 
[TOlwith X = —1 from a cyclic shift and Vervaat's theorem [ID]) that the law of 'i?"(c) under Pexc^ 
converges weakly to A^*-*^ = A^(- \ t = t). On the other hand, it can be checked that 4~"(n + 
l)^/nCatn increases in n towards its limit TT-^/^, SO that a3/22-2MCat^„t]Pexc*^(i^(i?"(c))) is 
dominated by K'e~'^^^ /\/t for some K' > 0. By dominated convergence, we get 

{2ay/^LMo i{Hor) ^ A I -^n'^^Xh) = AN{H) , 

by the fact that N{t G dt) = dt/V27rfi. Incorporating the spatial displacements is now easy 
from the convergence of the contour function, see for instance [23]. □ 

Finally, the measure associated with Q0.2 through the mapping ^'0,2 is the measure LMo,2 = 
2E^eMg=j LMm where 



LM„ = 5„,WN4dw;^*d^i;^*)DSi(d(c,z))#,(dr)DS,^._i(d(c^*,z^*))DS,^.(d(c^*,z^*)). 
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Its continuous counterpart is of course CLMo,2 = 2^^gjyj>3 CLMm, where 
and PNm takes the particularly simple form 



PNn,(du;^'du;^*) = / -^=PS^Q{dw'*dw'*) 
Jr+ v27rr 



We define the scaling operation by 

Proposition 4 For every continuous functional H satisfying 

\H\ < Kexp(-/3r[), 

for some K,P > 0. Then 

-\-lMo,2iH o CLMo,2(i?). 

The proof uses both the ingredients of the proofs of Theorem [5] and Proposition [3l but 
presents no new difficulty and is left as an exercise to the reader. 

5 Asymptotic enumeration in the non-planar case 

We want to derive Theorem [T] by a study of the generating function 

n>l 

which will be done by elementary singularity analysis. We will restrict our attention to the non- 
planar cases g > 1, as the planar case is already known. First note that we have the following 
alternative expressions 

qGQg {q,a^)GQg,i 

which is finally equal to LMg^i((r[/2-|-x(3))~^s'^'/^). The coefficients' asymptotic behavior is the 
same as that of Gq{s) := LMp^i(2r[~"^s'^"'^), so we may and will work with the latter. Now write 
s = e~^/", H = r[~^e~'^', and apply Theorem [3 We readily obtain that 



a—>oo ^ 

where C'g £ {0,oo) is defined by C'g = 2^3^ CLMg,i(Tf ^e"^')- Otherwise said, 

Go{s) ~ C'gil - s)-(5s-3)/2 (20) 
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as s — > 1 along [0, 1]. From this, Tauberian theorems allow to obtain 

n (-11 

which is weaker than (but consistent with) what we are aiming at. In order to apply standard 
transfer theorems of singularity analysis, we need to show that the asymptotics (j20p hold for 
s ^ 1 in a domain of the form 

^fl,<^ := {s e C : |s| < i?, s 7^ 1, 1 arg(s - 1)| > (/>} 

for some R> !,(/> G (0, 7r/2) and on which Gq has to be extendible analytically (here we take 
the determination of the argument in (— 7r,7r]). 

Lemma 12 Let a > and f{z) = X]„>o hnz'^ be a power series with non-negative coefficients 
satisfying hn ~ n'^'^ /T{a) as n ^ oo. Then f{z) ~ (1 — for z — > 1 inside a domain of the 
form Di ff, with (j) G (7r/2,7r). 

Proof. Let (1 — z)"" = X]„>o ^n-^") so that /i° ~ Fix e > 0, take N large enough so that 
— hn/h^\ < e. From the fact that /i^, /i^ > 0, we easily get 

N 



1(1 - zTfiz) - 1| < |i - zr Y.^K + hi)\zr + ef 



(] 

n=0 ^ 

and the first term goes to as z — > 1. As for the second, write 1 — z = pe^^ for \6\ < vr — (j). 
For p sufficiently close to 0, we get that |1 — z\/{l — \z\) < p/(l — (1 + — 2pcos(7r — (/)))^/^), 
which admits a finite positive limit as p — > 0, allowing us to bound the right-hand side of the 
last displayed expression by a constant (i;^-dependent) multiple of e. Hence the result. □ 

Now rewrite Go{s) as 

Go(s) = 2 /wNm(dr,dx) (g) DSrJ^s^'/A . 

We use the following trick: for a function H on M-^ we have 

where M^'^'° denotes unrooted maps and (m°,c*) is the map m° rooted at e*. Indeed, it holds 
that (m°,c) = (m°,e') as a rooted map for edges e, e' G E{xn) if and only if there exists an 
automorphism of m sending e to e'. Thus, the orbits of the natural action of Aut (m°) on E{m) 
correspond to the distinct possible roofings of m°. Moreover, if p • e = e, then p is the identity, 
which can be seen by reasoning in a step-by-step fashion from the edge e and going to neighbors. 
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Therefore, all the orbits have same cardinality |Aut(m°)|, and the claim follows. Applying it 
twice, once in each direction, allows us to get (with a straightforward definition for WNm°) 



Go{s) = 4 ^ |Aut(m°)ri /'wNn.°(dr,dx) (g) DS,, 



4 l^("^)r' / WNm(dr,dx) (g) DS,, 



By (Uni) we have DS,.(s^('=)) = ((1 - Vl - s'^)/sY, so that we get Go(s) = Gi o ^2(5), where 
Gi(z)=4 |^(m)r^ /"wNm(dr,dx)z^-6E(™)^% 



and 



s 

which can be extended analytically to a domain I?/?,^ with i? > 1 and G (0,7r/2). Moreover, 
elementary computations show that -D/j,^ is mapped to a domain Di g for some 9 £ (7r/2,7r). 

On the other hand, applying LemmaOto = 1 and pa = n = [(2a)^/^] gives that the n-th 
coefficient in Gi satisfies 

[z"']Gi{z) ~ 4 ( -Yn'^-^ y \E{m)r' [ = cy^-^ . 



Note that |£'(m)| = 2{6g — 3) for every term in the sum. From Lemma [T2l we deduce that 
Gi{z) ~ r{5g - 3)C^'(1 - ^)-(59-3) for z ^ 1 inside Di^g. Combined with the fact that ^2(5) = 
1 + 2\/l — s + 0(1 — s) as s ^ 1 in Dr^^, we finally obtain that Go has an analytic continuation 
in Dji^^, and is equivalent to 2~^^"'''^r(5g' — 3)G^'(1 — s)~(55'~3)/2 g _> j inside this domain. 
An elementary transfer theorem [181 Corollary VI. 1] allows us to conclude that [s"]Go(s) ~ 
G^n^*-^"^^/^, where the constant Gg is computed as 

C -i6py r(5g-3) ^ f AA(du) 

' W (6, - 3)r te) ^Ak ■ 



9,1 

Elementary (but tedious) computations give Gi = 1/24, consistently with [5]. □ 
Remarks. 

1. It is expected, and checked on numerous examples, see Gao |20j . that "reasonable" families 
of random maps with size n will have the universal property that their cardinality is equivalent 
to Ga^'n-^^^a)/^ for some G, a depending on the model. A nice feature of the above approach 
is that, let aside the prefactor 16(3/64)^, it is reasonable to think that the gamma terms and 
summation over ^ '^ot a particular feature of quadrangulations, and will appear for more 
general families of maps. However, it is extremely hard to determine the exact value of this 
sum as it involves a number of terms growing very quickly with g, and less and less tractable 
integrals. 
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2. In [To], the authors give an alternative expression for Cg, still in the form of a prefactor times 
a sum over trivalent maps with one face (and with an extra labeling of vertices). The expression 
of [To] is somewhat simpler than ours, as it does not involve integrals but only finite products of 
quantities that depend in an elementary way on the labeled map involved in the sum. It would 
be interesting to show directly that the two expressions are the same. 



6 Metrics on metric spaces 

In a metric space {X,d), we let Bxix,r) (or simply B{x,r) if there is no ambiguity) be the 
open ball centered at x with radius r. We let dn be the usual Hausdorff distance between closed 
subsets of X 

dH{C, C) = inf{e > : C C {C'f and C C C'] , 

where = {x ^ X : infyg^d(x,y) < e] is the e-neighborhood of A. We also let dp be the 
Prokhorov metric between Borel probability measures on {X, d) defined by 

dp{^i, /i') = inf{e > : /i(C) < ix'{C') + e for all closed C] . 

The topology associated with dp is that of weak convergence of probability measures on X. See 
|12^ Chapter 11] for references. 



6.1 Gromov-Hausdorff distance 

We first recall the definition and first properties of the Gromov-Hausdorff distance, taken from 
[221 El He]. If X = [X, d], X' = [X', d'\ G M, we write 



dGH(X,X') = ■mi 5H{(l>{X)A'{X')) , 

<p,<p 

where the infimum is taken over all isometries cj)' from X, X' into a common metric space 
(Z,5). 

An alternative useful description is given as follows. A correspondence between two sets 
X,X' is a subset 7^ C X x X' such that 7^(7^) = X and 7r'(7^) = X' , where 'k : X x X' ^ 
X^tt' : X X X' ^ X' are the canonical projections. We use xlZx' as an alternative notation 
for (x,x') G TZ^ and we let C{X,X') be the set of correspondences between X and X' . If 
(X, d), {X' ,d') are metric spaces, and 7^ is a correspondence between X and X' , the distortion 
of TZ is defined as 

disT^ = sup{|d(x, y) — d\x' , y')\ : xTZx' , yTZy'} ■ 

Then 

dGH(X,X') = ^ inf dis7^, (22) 

2 ■Rec{x,x') 

which of course does not depend on the particular representatives X,X' of X, X'. 

It is an easy exercise to check that disT^ = disT?. where TZ is the closure of 7^ in X x X' 
for the product topology, so that the infimum in (j22p could be taken over the set Cc(X, X') 
of compact correspondences (this set depends on the distances d, d' and not only on the sets 
X, X'). The following properties [HI [16], combined with the forthcoming Proposition |8l will be 
useful in proving Theorem [21 in Sect. 17.11 
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Proposition 5 (i) The function den is a metric on M. The space (M, cIgh) is separable and 
complete. 

(ii) A subset A C M is relatively compact if and only if the set {diam (X) : X G A} is 
bounded, and for every e > 0, there exists N > such that every X E A can be covered with N 
open balls of radius e. 

6.2 Gromov-Hausdorff-Prokhorov metric 

Let X = [X, d, /i], X' = [X', d', /x'] G We let 

dGHp(X,X') = inf {6h{4>{X)A'{X'))V Sp{(t>,f^,<py)) , 

where the infimum is taken over all isometries </>, (j)' from X, X' into a metric space {Z, S). 
Theorem 6 The function dcHP defines a distance on M„. 

This result is well-known, see for instance [441 p. 762], with the slight difference that the 
supremum between 6h and 5p in the definition of dcHP is replaced by a sum. In order to make 
this article self-contained, we will provide a full proof of Theorem [U to which the remaining part 
of this subsection is devoted. In order to handle the distance dcHP, we develop an alternative 
definition with the help of correspondences and couplings. This is inspired in part from [221 121j . 

If n, n' are Borel probability measures on (X, d) and (X', d'), we say that a Borel measure v 
on the product space X x X' \s a. coupling between /i and fi' if vr^i/ = /i and vr^z^ = /i'. We let 
Ai{fi,fj,') be the set of couplings between /z and fi'. It is never empty as it contains the product 
measure and is closed (and even compact since X x X' is compact) for the weak topology. 

A useful feature of couplings is the following. Since all spaces that are considered are Polish, 
it is known that we can write a coupling z/ G M.{n, fi') in the form n(dx)Qi,{x,dy), where Qi, 
is a Markov kernel from X to X'. In particular, if ^ui, . . . , ;Ufc are probability measures on the 
compact spaces Xi, . . . , Xf^ and i^i,! < i < k — 1 are couplings between //j and /ij+i respectively, 
we can define a Markov chain {Ai, . . . ,Ak) such that Ai is a random variable with distribution fii 
on Xi, with initial distribution ni and transition kernel Q^^ at step i. The joint law of (^i,^^) 
is then denoted by i^i . . . i^^ and is a coupling between /xi and /z^. 

Proposition 6 Let X = [X,d,fi\,X' = [X',d',n'] G M„. Then 

Proof. Assume that dGHp(X, X') < e and let (j), (j)' be isometric embeddings of X, X' into (Z, 5) 
such that cj)'{X')) < e and 5p{(j)^fi, 4>'^^') < e. Arguing as in [21 Remark 7.3.12], we may 

and will assume that (Z, 5) is separable. Then the set TZ = {{x,x') G XxX' : 5{(f){x), (p'{x')) < e} 
is an element of Cc{X,X'). To evaluate its distortion, let xTZx' and yTZy', and write 

d{x,y) - d'{x',y') = 5{<P{x),^{y)) - 6{cP' {x'), <p' {y')) < 6{cP{x), cP' {x')) + 5{<P' {y'), cPiy)) < 2e , 

and the lower bound is obtained in a similar way. Thus disT?. < 2e. 
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Now, by a characterization of the Prokhorov distance (see \12\ Corohary 11.6.4]) on separable 
metric spaces, we can find a coupHng D E A4{(j)^fj,, 4>'^ii'), which is a measure on Z x Z, such that 

i'{{{x,y) : S{x,y) > e}) < e . 

We let 1^1 S M.{fj,, (p^fi) be the image measure of /i by j; i-^ {x,cj){x)) and 1^2 G -^((A*^'! /^') be 
the image measure of fi' by x' 1— > {(j)'{x'),x'). Then the coupling u = vivu2 obviously satisfies 
v{TZ'^) = u{{{x^y) : 5{x,y) > e}) < e (here TIP is the complementary set of TZ). 

Conversely, assume given TZ € Cc{X,X') and u G (//,//') as in the right-hand side of (f^Sj) 
for the value e. As in [9] we endow XL\ X' with a metric 5 that is compatible with d, d', i.e. its 
restrictions to X and X' coincide with d,d' , and such that 5{x,x') = e whenever xTZx' . To do 
this, we set for every a; € X, x' G X' 

6{x, x) = \ni{d{x, y) + e + d'{y', x') : yUy'] . 

The only non-trivial point to show is the triangular inequality, and we prove only the hardest 
case, the others being symmetric or left as an exercise. Let x,y £ X and x' € X', and consider 
zTZz',tTlt'. Since by definition \d{z,t) — d'{z',t')\ < 2e we estimate 

{d{x, z)+£ + d'{z\ x')) + {d{y, t)+£ + d'{t', x')) 

> {d{x, z) + d{z, t) + d{y, t)) + {d'{z', x') + d'{t', x') - d'{z', t')) 

> d{x,y)+0, 

which is obtained by two applications of the triangular inequality in X and X' . Taking the 
infimum over (z, z'), {t, t') allows us to conclude that d{x, y) = d{x, y) < 6{x, x') + 5{x' , y). 

Now u induces a coupling v on the disjoint union X lA X' between i*^ and i'^/i', where 
are the canonical injections X ^ X L\ X' and X' XU X' . Now obviously for every a > 1 

i>({(x,x') : 5{x,x') > ae}) < z^(7^^) < e, 

so that 6p{i*fJ-,i'^fJ.') < £• Hence the result. □ 

We are now ready to prove that dcHP is indeed a distance. Symmetry is obvious, as well as 
dGHp(X,X) =0. 

Separation. Let X = [X,d,ij],X' = [X',d',fi'] G M„, with dGHp(X,X') = 0. Then for every 
n > 1 we can find Tin G Cc{X,X') and E M{fi, n') with dis7^„ < and z^n(7^^) < 2"'". 

Endow the product space X x X' with the distance S{{x, x'), {y, y')) = max((i(x, y),d'{x' , y')). 
Then {X x X',5) is a compact metric space, so that the set of compact subsets of X x X' , 
endowed with the Hausdorff distance 6h, is compact [9l Theorem 7.3.8]. Therefore, {TZn,n > 1) 
converges up to extraction to some compact TZ. It is immediate to check that the latter is an 
element of Cc{X, X'), with disT?. = 0, and this implies that TZ is the graph of a bijective isometry 
h:X^X'. 

On the other hand, since A4(//, fi') is compact for the weak topology, the couplings (f^, n > 1) 
converge weakly, up to re-extracting, to some G A^(/.f,/i'). Since TZn converges to TZ for the 
distance 6h, it holds that for every e > 0, TZn '^TZ = {{x,x') : inf(j^ y-jg-;^ (5((x, x'), {y,y')) < e} 
for every n large enough. For such n, we have Vri{TZ^) ^ ^niTZn) > 1 — 2~", and by well-known 
properties of weak convergence of probability measures, since TZ is closed, this entails that 
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1^(7^^) > limsupt'„(7^'^) = 1. Letting e ^ shows that = 1, and we conclude that is 

supported on the graph of the isometry h. Thus, for every measurable g : X' ^ M_|_, 

/ g{h{x))fj,{dx) = / g{h{x))i>{dx,dx') = / g{x')iy{dx,dx') = / g{x')fi'{dx'), 
Jx JXxX' JXxX' Jx' 

so fi' = h^:^. Since X' = h{X), we obtain X = X'. 

Triangular inequality. Let X = [X,d,fi],X' = [X' , d' , fi'],X" = [X",d",ii"] G M„, and let us 
assume dGHp(X, X') < ei and dGHp(X',X") < £2- We wish to show that dGHp(X, X") < £1 +£2- 
By Proposition El we can find TZi G Cc(X,X'),7^2 G Cc{X',X") and 1^1 G M^fi, fi'),i^2 G 
M{n',fj,") such that dis7^i < 2ei and ^'^(7^f) < ei,i G {1,2}. 

Let 1^3 = 1^11^2 and define a correspondence TZ^ G C(X, X") by xTZ^x" if and only if there 
exists x' G X' such that xTZix' and x'TZ2x" , that is 7^3 is the image by the (continuous) canonical 
projection X x X' x X" X x X' of the compact set {TZi x X") n (X x 7^-2), and hence is 
compact. To evaluate the distortion, let xTZix\x'TZ2x" and yTZiy' ,y'TZ2y" and write 

\d{x,y) - d"{x",y")\ < \d{x,y) - d'{x',y')\ + \d'{x',y') - d"{x",y")\ < 2(£i +£2) , 

so that dis7^3 < 2(£i +£2). 

Also note that l(7^3)c(2;, x") < lniix,x') + ln^^{x',x") for every {x,x',x") £ X x X' x X" . 
Indeed, if {x,x") ^ 7^3, then for all x' G X', either {x,x') ^ TZi or {x',x") ^ 7^2- We deduce 

1^3(^3) < / fiidx)Qi{x,dx')Q2{x',dx")lTii{x,x') 

JxxX'xX" 

+ [ ii{dx)Qi{x,dx')Q2ix',dx")lndx',x") 

JxxX'xX" 

= (7^^ ) + 1/2 (7^2 )< £1 + ^2 • 

This is the wanted result, concluding the proof of Theorem [6l □ 
6.3 Relation with other metrics 

It turns out that the topology determined by dcHP has already been considered by [IHllIT], where 
a different distance was introduced. Note that in [T7j, the authors restricted their attention to 
compact spaces which are M-trees, but all of their results up to Theorem 2.5 in this reference 
can be rewritten for the space without any change. 

An £-isometry from X to X' is a mapping / : X — > X' such that {f{X)Y = X' (we say that 
f{X) is an £-net in X') and for every x,y £ X, 

\d{x,y)-d'ifix),f{y))\<s. 

We say that X = [X,d, ij],X' = [X',d',i2'] G M„ are £-close if there exist two measurable s- 
isometries / : X — > X' and 5 : X' — > X such that d'p{f^fi, fj,') V dp{g^fj.', n) < e. Obviously, being 
£-close is indeed an isometry class property and does not depend on the choice of representatives. 

We define Aghp(X,X') as the infimal £ such that X, X' are £-close. This does not define a 
true distance, as the triangle inequality is satisfied only up to a factor of 2. This quasi-distance 
is however separated and thus sufficient to define a separated topology on Mw It is also fully 
discussed in [TTj, building on metrization theorems for uniform spaces how to build a true 
distance using AghPj which is yet different from dGHP- 
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Proposition 7 It holds that 3 ^Aghp < ^gup < 2Aghp- 

Proof. Assume that dGHp(X, X') < e, and let 7^, be a correspondence and a coupling between 
X and X' such that disT^ < 2e and > 1 — £• We build a measurable 3e-isometry f : X ^ X' 
out of TZ, by taking {xi, . . . ,xj\f} a minimal e-net in X , choosing for every i some such that 
XiTZx[, and finally defining /(-Bj) = {x^} where i?i = B{xi,e) \ \J^^j^^_j^ B{xj,e). Then the 
image measure P of by (x,x') i— > {f{x),x') is a coupling between /^.^ and Moreover, from 
the way / is constructed, if d'{f{x),x') > 3e then {x,x') ^ TZ, so that z>{(x',y') : d'{x',y') > 
3e} < v{TZ^) < e. Thus d'p{f^i_L, fi') < 3e, and by a symmetrical argument we conclude that X 
and X' are 3e-close. 

Conversely assume given f : X ^ X' an e-isometry with d'p{fi:fi, fi') < e. Define a corre- 
spondence between X and X' by xT^x' if d'{x' , /(x)) < e. Its distortion is estimated as follows: 

\d{x,y)-d'{x',y')\ 

< \dix,y) - d(/(x),/(y))| + \d{f{x),f{y)) - d{f{x),y')\ + \d{f{x),y') - d{x',y')\ , 

which is less than 4e whenever xTZx', ylZy' . We take the closure of TZ, which has same distortion, 
and still call it TZ. Take a coupling i> between f^fi and /i' with ^{(x', y') : d'{x' , y') > e] < e and 
let vi be the image measure of /U by x i-^ (x,/(x)). Then = uiu satisfies v(TZ) > 1 — e. We 
conclude that dGHp(X, X') < 2e. □ 
As a consequence of this and \n\ Lemma 2.3], we obtain that the metric properties of Evans 
and Winter's distance are the same as those of dGHP- In particular: 

Proposition 8 (i) The metric space (M„,dGHp) is separable and complete. 

(ii) A subset A of M„ is relatively compact if and only if the subset { [X, d] : [X, d, /.i] € A} 
is relatively compact in (M, dGn)- 



6.4 Marked topologies 

In this section we introduce a marked variant of the Gromov-Hausdorff topology. 

A fc-marked metric space is a triple (X, d, (Ci, . . . , Cfc)) where {X,d) is a compact met- 
ric space and Ci,...,Ck are compact subsets of X. The isometry classes of marked spaces 
[X, d, (Ci, . . . , Ck)] are defined in the obvious way, and their set is denoted by M^. We let, with 
obvious notations, 

d|H(X,X') = inf (6h{HX),cI>'{X')) V max <A'(CD)) , 

4>i4> \ i<i<k J 

the infimum being taken over all isometric embeddings ^ of X, X' into a common metric space 

The following statement is a straightforward generalization of the classical results of Propo- 
sition [5l and its proof is left as an exercise to the reader. 

Proposition 9 (i) The function dgjj is a distance on which is alternatively described by 

d^H(X,X') = ^infdis7^, 
the infimum being taken over correspondences TZ G Cc{X, X') such that TZr\{Ci x C^') G C(Ci, C'^). 
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(ii) The metric space (M^,,dQjj) is separable and complete. 

(iii) A subset A C is relatively compact if and only if {[X, d] : [X, d, (Ci, . . . , Ck)] G A} 
is relatively compact in (M, cIgh)- 

When Ci,...,Cfc are singletons, we write for notational convenience {X,d,x.) instead of 
{X, d, ({xi}, . . . , {xfc})). The reader will easily be convinced that the subspace of such marked 
spaces is closed in M^. 

6.5 Randomly marked spaces 

We now introduce random marks on metric spaces. If X = [X, d, fj] € M^, we can consider an 
i.i.d. sequence (C/i, . . . , Uk) which is /i-distributed in X, and we want to consider the isometry 
class of [X,d, (Ui, . . . , C/fc)] as a random variable in M^. We also wish to consider randomly 
marked versions of spaces that are themselves random. To do this implies that we check a 
couple of measurability issues. We introduce for each k a kernel from to M^, defined as 
follows: 

mfc(X,A)= / ^®^(dx)lA([X,d,x]), X=[X,d,fi\eM^,Ael3{Ml). (24) 

Note that the integral makes sense as the function X'^ — > which maps x to [X, d, x] is 
obviously continuous, so that x i-^ 1a([-'^, d, x]) is measurable by composition. Also notice 
that it is independent of the particular choice of the representative {X, d, /i) G X, and hence is 
unambiguous. 



Lemma 13 The formula \24^ defines a Markov kernel. 

Proof. It is obvious from the definition that for every X, the function A nifc(X, A) defines a 
measure on M^. 

Now let A be closed in M^. We show X mfc(X, A) is upper-semicontinuous. Taking e > 0, 
for any X' G with dGHp(X, X') < e, and with obvious notations, we can find TZ € Cc{X,X') 
and V G M.{^,ii') such that > 1 — e and disT^ < 2e. We write 



f ^®^(dx)(l-lA.([X,d,x]))= I I.^'=(dx,dx')(l-lA^([X,d,x])) 
J X*' J{XxX')'= 

< i-(^i-ef+f i/®^Xdx,dx')l{,,7^x.M<i<fc}(l-lA^([^,t^,x])) 

< l-{l-ef+f d(^'r'(x')(l-lA([^',d',x'])) , 

J{X')k 



where we used the fact that j^®'^({(x, x') : Xi7lx[, I < i < k}) > (1 — e)*^ and that whenever 
the latter event occurs, it holds that dQjj([X, d, x], [X',(i',x']) < e. Since A^ decreases to A as 
e J, we see that the left-hand side can be made arbitrarily close to 1 — nifc(X, A) provided e is 
sufficiently small. Therefore 

limsupinfc(X', A) <iiifc(X, A), 

as claimed. We conclude by a monotone class argument that mfc(-, A) is measurable for all Borel 
A C M^. □ 
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A random variable with law mfc(X, •) can be considered as the rigorous definition of the space 
X randomly marked with k independent /u-chosen points. If M(dX) is a distribution on M^, a 
random variable with law M(dX)nifc(X, •) is interpreted as a weighted space chosen according to 
M, and then marked randomly using the weight. 

A key property of randomly marked spaces is the following. If X G M^, we let X'^ be the 
canonical random variable with law nifc(X, •). A similar notation holds if X is a random variable 
with values in M„. 

Proposition 10 Let (Kn,n- ^ 1) be random variables in M„ converging to X in distribution. 
Then the randomly marked spaces X^ converge in distribution towards X'^, for the topology of 
M*. 

Proof. Using the Skorokhod representation theorem, we may assume that the probability space 
(0, P) on which the random variables are defined is such that the convergence X„ X holds 
a.s. Now, we have from the proof of Lemma [T^ that mfc(-. A) is upper-semicontinuous for every 
closed A C M^,. By the reversed Fatou lemma this gives 

limsup^[mfc(X„, A)] < E[mk{X,A)] , 

n—foo 

which says that lim sup ^(X^ G A) < P(X'^ G A), and the conclusion follows from a well-known 
characterization of weak convergence. □ 

6.6 Three lemmas 

We end this section with key lemmas to the proof of Theorem [3l 
6.6.1 Diffuse measures 

Lemma 14 Let d„, n > 1 converge m towards [X,d,fj]. Assume that 

limliminf / fin{<ix) ^n{B {x , e)) = . 

Then ^ is a diffuse measure. 

Proof. From the proof of Lemma fTSl we know that if A is an open set in M2, then X 1— > m2(X, A) 
is lower-semicontinuous from M„ to Mg. Now write 

/ Hn{dx)fIniB{x,e)) = n'^'^{dx,dy)l[d„(x,y)<e} =M^n,A.) , 

Jx„ Jxl 
where A = {X' = [X\d', {x' ,y')] G : d'{x',y') < e} is open. Therefore, it holds that 

liminf / /i„(dx)/i„(5(x, e)) > / fi{dx) ^j,{B {x , e)) . 
Jx„ Jx 

Letting e — > and using the hypothesis gives the result. □ 
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6.6.2 Support 

Lemma 15 Let d„, /i„,], n > 1 be a sequence converging to [X,d,fi] in M^. Let Yn X, 
be En-dense in Xn, for some — > (i.e. every point of X^ is ctt distance at most of a point 
ofYn). Assume that for every e > 0, 



n 



limsup inf fj.n{B{x,e)) > 0. 



Then supp fj, = X . 



Proof. Take x & X and e > 0. We must show that fi{B{x,£)) > 0. To this end take compact 
correspondences TZn between Xn,X with distortion going to 0, and couphngs between 
and /i with UniJ^n) — > 1- Let Xn be such that Xn'R-nX and x'^ £ be at distance < e„ from Xn- 
Then 

fl[B{x,e))= I tJ-{dy)l{d{x,y)<e} > / l'n{<^ynAy)^{d{x,y)<e} ■ 

Jx JXnXX 

Now take n large enough so that disT^.^ < e/2 — e„. Then dn{x'n,yn) < e/2 implies d{x,y) < e 
as soon as ynT^nV- Therefore the last integral is greater than or equal to 



l^n{dyn,dy)l{d„{x'„,y„)<e/2} > / M'^Vn, dy)l{d„(^^>^,y„)<e/2} - '^n{T^n) 

= ^„(B(4,e/2))-z.„(K)- 

This is bounded below by inf^^gy^ Hn{B{x, e/2)) — f„('7^^), which has positive limsup as n — > oo. 
Hence the result. □ 



6.6.3 D-median points in geodesic metric spaces 

Since compact geodesic metric spaces form a closed subset PM of (M, cIgh) ([HI Thorem 7.5.1]), 
it is immediate that the set PM„ of spaces [X, d, G for which [X, d] is a geodesic metric 
space is closed as well, and a similar statement holds for /c-pointed geodesic spaces, whose set is 
denoted by PM;^. 

Let {X,d) be a compact geodesic metric space, and let x,y G X and D £ {—d{x,y),d{x,y)). 
A point z € X is called a D-median point between x and y if d{x, z) — d{y, z) = D. Note that 
the set med^j^(X) of such points is compact and separates the space X into two open subsets, 
one containing x and the other y, since every path c : [0, 1] — > X from x to y is such that 
d{c{t),x) — d{c{t),y) ranges from —d{x,y) to d{x,y) continuously. We also say that z is on a 
5-quasi geodesic from x to y if d{x, z) + d{y, z) < d{x, y) + 6. We let qg^^y{X) be the set of such 
points, and medqgf/(X) = medfj^(X) n qg^y(X). 

Proposition 11 Let [Xn, dn, (x„, G PMg be converging to [X, d, {xi,X2)], and let also Dn G 
{—dn{xn,yn), dn{xn,yn)) be Converging to D £ {—d{x,y),d{x,y)). We assume that for every 
e > 0, we can find 5 > such that, 

limsup diam (medqgf;^j;f^(X„) ) < e. 

n— »oo ^ ^ 

Then there exists a unique element of med^y{X) which is on a geodesic path from x to y, i.e. 
m.edqg^y^ {X) is a singleton. 
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Proof. Assume that there are two possibly different points z,z' € med^y{X) that belong to 
some geodesic path between x and y. Choose arbitrarily such a path 7 = 'y^y going through z 
and let e > 0. Let e > (5 > be as in the hypothesis for this value of e, and let = io < < 
. . . <tK = d{x, y) be such that d{u^ , uf'^^) < 5/7 where = 'j{tk)- Let TZn be a correspondence 
between X„ and X such that XnTZnX and yrJ^nV-, and dis7^„ — > as n — > 00. We may assume 
that n is large enough so that dis7^„ < 5/7. Choose arbitrarily v/^TZnU^ and let 7„ be the 
concatenation of arbitrarily chosen geodesic paths between v/^ and v!^^ for Q < k < K — 1. 
Take a point v on 7„. Then for some k, v lies on a geodesic from vl^ to lijj,"''^. Then 

< (i„(n^+\u^) + (i„(n^,Xn) 

< d{u''+^,u'') +d{u'',x) + 2disnn 

< x) + 2dis Tin + 5/7 < d{u'' ,x) + 35/7, 

so that by a symmetric inequality involving yn and by summing, 

dn{v, Xn) + dn(?^, Un) < d{x, y) + 65/7 < dn{Xn,yn) + ^ • 

Therefore the image of 7^ is included in qg^^y^{Xn)- Let Zn be any element of med^J^y^{Xn) 
that belongs to 7^, and let £ such that ZnTZnZn- 

Now doing the same with a geodesic path 7' from x to y passing through z' , we can construct 
elements z'JZnz'^ such that z'^ G medf;'j^^^(X„) n qgi^^y^{Xn). Hence dn(^n,4) < e for large n 
by hypothesis, and thus d{zn, z'^) < e + dis7^„. 

Letting n go to infinity, possibly through a subsequence, we obtain that z^ — > z{e) and 
z'n z'{e), where necessarily z{e),z'{e) G med^y{X) by passing to the limit in 

\d{zn,x) - d{zn,y) - Dn\ < 2dis 7^„ , 

and similarly for z'^. By a similar argument, since z„ G 7n we have z{e) G 7*^, the 5-neighborhood 
of the image of 7. Since 7"^ decreases to 7 as 5 | and med^j^(X) is compact, letting e (and 
5) decrease to 0, z{e) converges to the unique element of 7 n med^y{X), which is z. Similarly, 
can be made as close to z' as wanted. But since d{zn,z'^) < e + dis7?.„ we get that 
d{z{e), z'{£)) < £, and this entails z = z'. □ 

7 Metric aspects of random quadr angulations 

We now embark in the proof of Theorems E] and O We start with a preliminary remark. For 
technical reasons, the set PM and its weighted and marked analogs are somewhat easier to work 
with for our purposes. Of course, the main space of interest, namely Xq for a quadrangulation 
q, is not a geodesic metric space. However, its dGH-distance to PM is at most 1. Indeed, take 
copies of the unit segment (with the usual metric) {x*^ : < 2; < 1} indexed by an orientation 
-^1/2(^1) of the edges of q, and identify the points 0^ for which e" coincide, and the points 1*^ for 
which coincide. The quotient metric graph [9l Chapter 3.2.2] is a geodesic metric space whose 
restriction to the points 0^,e G £'i/2(q) is isometric to {V{q),dq), and obviously its Gromov- 
Hausdorff distance to the latter is less than 1. Thus any limit in distribution of a"^/^Xq must be 
at distance from PM with probability 1, hence belong to PM. In passing this proves property 
1. in Theorem [3l For these reasons, we work with the metric graph constructed above, which 
we still call Xq, and restrict our attention to geodesic metric spaces from now on. 
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7.1 Proof of Theorem m 

For X G M and e > 0, we let AA(X, e) be the minimal number of open balls of radius < e that 
are needed to cover X. For any sequence N = {N^, A; > 1) of positive integers, we introduce the 
set 

IC^ = {X™ = [X,d,fi] G PMw :J\f{[X,d],2^'') < A^fc Vfc > 1} . 

By Propositions E] and and ^ Exercise 7.4.14], the sets ICjv are (relatively) compact. By 
Prokhorov's theorem, it is thus sufficient to show that for every e > 0, there is a choice of A^, A 
for which 

Qf ({a-i/'^X^ G Kn, < a- < a}) > 1 - e . 

To do so, we use the results of Sections [SHU in the case k = 1. We also separate the case g > 1 
and 5 = 0, starting with the former. Consider (q, x) G Qg,i and let (m, [1]) = ^'g^i(q, x) and 
(m, i,t*) = ^g^i{ci,x). Recall also the definition ([TT]) of /, and write 

uj{h,e) = sup \h{x) — h{y)\ 

x,y£j,\x—y\<e 

for the modulus of continuity of a function h defined on an interval J. 
Lemma 16 For R > 0, 

{AA(Xq, R)>N}C [2u; (l, [^j ) > r} (25) 

Proof. Since m has only one face of degree T[, all the half-edges of m are of the form e{j) = 
Vm(c*(0)) for some integer j G [0, r[) with the notations of Sect. 13. 3i Split the lifetime interval 
[0,Ti] of / into N disjoint parts Ii,...,In of lengths at most \ti/N~\. By Lemma E] and the 
definition (jlOj) of I, if uj{l, \ti/N~\) < R/2, then all the vertices of m of the form e{j)~ for j G li, 
and considered as vertices of q, must be within dq-distance at most R, for each 1 < i < A^. 
Since all vertices of q are contained in this way except the distinguished vertex x, adding yet 
another ball to cover x yields (|25|). □ 

We now consider a marked version of Q^^^ as the probability on Qg,i defined by 

Q„i(V-ie-^^<id(q,a;)) 

Using Section m its image under ^'^.i can be written 

LM,,i ((V2 + x{9)r' exp (-/3V2) d(m, [, t* 
LMg 1 •= ^ 



LM,,i (^(rr/2 + x(5))" exp(-/3ri/2: 

Now, Theorem (and elementary bounds to get rid of the term xid)) entail that the law of 
e''(m, l,t*) under IM^P converges weakly towards 

CLm(^) ■= CLM,,i jr^'expi-Pn) d(m,[,r)) 
^'^ ■ CLM,,i(rfiexp(-/3r[)) 
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By the continuous mapping theorem, the distribution of (^/^"(Z^), e G E[m)) under LM^^l""^ jointly 

converge towards the law of the processes {l\ e e E{m)) ([H]) under CLM^'^^', and thus the law of 

under IM^i"^ converges towards the law of I under CLM^'^j^. Combined with the fact 
from Lemma [T6l that 

{AA(a-i/^Xq, 6)>n]c {2u; (v^"(/), (la)-' [^] ) > 6/,] , 



and by usual relative compactness criteria derived from the Ascoli-Arzela theorem [5l Theorem 
8.2], we see that for every k > 1, there exists A^^ such that for every a, Qg^^'^\M{a~^^^yiq, 2~^) > 

Nk) < e2-^. This entails Qf /"V"^^^Xq G K^) > 1 - e for this choice of = (A^fc, k > 1). We 
see as well that the random variables V^/a are tight in (0, oo) as they converge in distribution 
towards the law of T[ under CLM^^l. □ 

The case g = is similar but simpler, the role of CLMo,i being performed by the Ito measure 
of the Brownian snake, and taking into account the factor appearing in the definition of Qq''^^ , 
we must replace the limit in law of by the law of z under 

CLMo,i(r exp(-/3r) d(c, z)) 
CLMo,i(rexp(-/3r)) ' 

One will notice that the denominator would diverge if it were not for the corrective term r in 
front of the exponential, since CLMo,i(r G dt) = 4{27rt^)-^/^dt. The statement of Lemma [H] 
remains true without change, and we conclude in the same way. 

7.2 Proof of Theorem H 

Property 1. of this theorem has been discussed at the beginning of this section. Let us set things 
up for Property 2. We give the proofs for g > 1, the case g = being similar and easier. 

7.2.1 DifFuseness of fj, 

By Theorem[21 there exists an increasing sequence (a^, n > 1) along which the distribution under 

converges to a limit <^g^\ From now on, we will always suppose 
that a is taken along this sequence. From this and the discussion above (|26p . we obtain that 
the laws of (a~^/^Xq , a~^lq, ^'^(5'g_i(q, x))) under Qg^^"'*(dq)/iq(dx) form a relatively compact 
family of probability distributions on PM^ x (0, oo) x Cg i. Up to a further extraction, and by 
a use of Skorokhod's embedding theorem, we may and will assume that we are working on a 
probability space {i},J-,P) under which are defined random variables (qa,a^a) with respective 
laws Qlf^'^\dq)fiq{dx) such that 

1. letting = [Xa,da,iJ.a\ ■■= X^^, we have a"^/^X™ ^ X"^ = [X,d,fj], a random variable 
with values in PM^ (and we let X^ = [Xa, da]), 

2. letting (nia, [1^]) = ^'g,i(qa,Xa) and {ma,ia,t*a) = $g,i(qa,Xa), then ^"(ma,[a,il) ^ 
(m, l,t*), a random variable with law CLM^'^^^ 
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these convergences holding almost-surely as a ^ oo (along the proper sequence). We will adopt 
the notation la = {w^, c^, z^, z G E{m.a))^ = -^a + ° ^ai l^t be the duration of c^, and let 
la be the process defined as in ([TT]) or (fTO|) . Note that {2a)~^T\^ and a~^Vq^ both converge to 
r[ as a — > oo, so the random variable (X™,T[) has law ■5^g^\ By the previous convergences, the 
processes ip^'i^a) converge to so that ^"'{la) converges a.s. as a — > cx) towards the concatenated 
process I associated with (m, l,t*). 

Lemma 17 We have 

Ha{dx)fia{B{x,e)) 



J a 



2 



< (^) / dnd7;l||^a(;^)(L2a«J/(2a))-V-aa)(L2a^;J/(2a))|<eA} +4t^/- 

Proof. With each integer time i corresponds the exploration of a corner e{i) = V'ma (^* 
vcia, and each corner of ma is explored in this way because m^ has only one face. Considering 
the vertices e(i)~ as vertices of qa, we have da(e(z)~, e(j)~) > \la{i) — la{j)\ by LemmaEl From 
this, it follows that 



where the last term amounts from the fact that the selected vertex Xa is the only vertex of 
not corresponding to a vertex of m^, and a rough upper bound. We conclude by dividing by 
(t(^/2 + x(5))^ = ^„ ^iid a linear change of variables. □ 

Lemma 18 Almost-surely, the level sets {t G [0, r() : l{t) = x} have zero Lebesgue measure for 
every x G M. 

Proof. From the construction of / as the concatenation of the paths and given m, it suffices 
to show the same property for the processes = z^ -\- o cf of (I12p . Under CLM^^i, recall 
that conditionally on {w^, c^, t G E{m)), < t < r^) is a Gaussian process with covariance 

Gov {z^{s), z^{t)) = inisAt<u<svt c^{u), by definitions of snakes. Therefore, under this conditioned 
measure, the laws of pairs {l^{s),l^{t)) have densities, and consequently, for fixed s,t, 

ClMg^i{{l'{s) = l'{t),T' > s V t}) = , 

and the same is true for CLM^^l instead of CLM^.i. By Fubini's theorem, and with a slight abuse 
of notation, this implies that for any T > 0, ii U,V are uniform, independent random variables 
on [0, T] independent of (m, (, t*) then it holds that 

CLMf 1{{1'{U) =1'{V),t' > T}) =0. 

But obviously, if some level set {/^ = x} had positive Lebesgue measure with positive probability, 
then the latter probability would be positive for some T > 0. □ 

We now conclude that n is diffuse. Indeed, by the convergence of '4'"' {la) to / and of T[^/(2a) 
to T[ G (0, oo), we obtain from Lemma [T71 and dominated convergence that 



liminf / /ia(dx)^a(S(x, e)) < ^ / dudt;l||;(„ 



[0,r[)2 
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The last integral decreases to diid'yl{/(„)=/(„)}. as e — > 0, and this is by Lemma [T8l 

Lemma O allows us to conclude. 

7.2.2 Support of ^ 

The setting is the same as in the previous section. For e G £'(01^) we let x{e) = if 
dmai^^ J ^(n^^^)) > '^ma(e~, ^(m^^)) and x{e) = e~ otherwise (intuitively, it selects the vertex 
incident to e that is the furthest away from m-^). For u £ [0, T[^) we let 



[u\ ifx(e(H))=e(LnJ)- 
\u] ifx(e(H))=e(LnJ)- 



Lemma 19 The image measure o/2~^dul[o_T-[ )('") by u ^ e{{u)a)~ is the measure on V{ma) = 
V{qa)\{xa} assigning mass 1 on every vertex of V {ma)\V (ma) , where ma is seen as a subgraph 
of rWa, and assigning mass deg^^(x)/2 to every vertex x G V{ma), the latter denoting the degree 
in the graph ma- 

Proof. Assume x is an element of V{'nia) \ ^(m^^). Then e{{u)a)~ = x if and only if i = [nj is 
either the first integer such that e(i)^ = 2: or the last integer such that e(i)~ = x, corresponding 
to the two visits of the parent edge of x in the tree rooted on ^(m-^) containing x. The total 
mass of such u's is thus 1/2 + 1/2 = 1. 

Next, assume that x G V{in^^) \ V{ma)- This time e{{u)a)~ = x if [u\ is one of the two 
integers i such that e{i) £ E(m^^) and e{i)~ = x. 

Finally, if x G V{ma) then e{{u)a)~ = x if \u\ is an integer i such that e{i) G -E(m^^) and 
e{i)~ = X. There are deg^^(x) such integers. □ 

The measure appearing in Lemma[l9]is of total mass T(^/2, and we let fla be the probability 
measure obtained by renormalizing it by this quantity. Obviously, the Prokhorov distance in 
a~^/^Xa between fia and jla vanishes as a — > 00, so that [Xa,a-^/'^d a , jla] converges to the same 
limit [X,d,fj] as [Xa,a~'^/^da, fJ^a]- 

Lemma 20 It holds that 

2a /■^'a/(2a) 

X&V{cia)\{Xa} SG[0,T[^/(2a)) Ti^ Jq 

Proof. Let x be a vertex in V{c{a) \ {xa}-, and let ix be an integer such that e{ix)~ = x. If j is 
another integer time, then da{x, e{j)~) < laij)+Liix) - 2 minjAi^<u<jvi,; L{u) < '^^{lai \j - ix\) 
by Lemma [3l Thus 

(tiJ2)/2, ({y G y(qa) : 4(x,y) < a^'^e}) = ^ ^'"'^ d^x l{d„(x, e((n)J-) < a^'^e} 

> If ''^ dul{oj{la,\{u)a-ix\) <a'^''e/2}, 
^ Jo 

where we used Lemma [19] in the first step. The result is obtained by taking the infimum over x 
and by simple rearrangements. □ 
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Since '<p"'{la) converges uniformly to the continuous limit Z, there exists rj such that for every 
(large enough) a, ij{'tp°'{la),i]) < e/(4^). From this, we readily obtain from the last lemma that 

limsup inf /ia(5^-i/4^ (x,e)) > 7?/2 > . 

Since V{cia) \ {xa} is a~^/^-dense in X^, one concludes by Lemma \TE\ that fi is of full support. 
7.2.3 Uniqueness of geodesies 

We now specialize the results of Sect. [5] to the case k = 2, starting with some general remarks. 
A delay D = [di, ^2] G ^?(q, (xi, X2)) is identified with the quantity di — d2, which is an integer 
in the interval (— iiq(xi, X2), (iq(a;i, X2)) such that D + dq{xi,X2) is even. 

Let (q, {xi,X2), D) £ Qg,2 and (m, [1]) = ^p_2(qj (xi, X2), D). Consider any chain in q from 
xi to X2 with length dq{xi, X2)- There must be a (unique) vertex x in this chain that belongs 
to Vi n V2 (where Vi = Vi{ci,x, D)), i.e. for which dq{xi,x) — dq{x2,x) = D. We claim that 
l(x) = min\/jny2 Ky)- Indeed, we know from Theorem H] that dq{y,xi) = l{y) — l(xi) and 
dqiy, X2) = l(y) — 1(^2) for y £ ViH V2, and since dq{x, xi) + dq{x, X2) = dq(xi, X2) since x is on 
a geodesic, there must be equality in the triangle inequality 2l{y) — l(xi) — l(x2) > 'iq(xi,X2), 
and l(x) is indeed minimal. In passing, we have obtained the formula 

(iq(xi,X2) = 2 mi_n 1(e) - l(xi) - I(x2) 

ee/i:ee/2 

= 2 min 1(e) — minl(e) — minl(e) — 2 

ee/i:eG/2 eS/i ee/2 

= 2 min_ 1(e) — min 1(e) — min 1(e) — 2 , (27) 

where is the face of m incident to its root e^,, and is the other face. Indeed, with a vertex 
in Vi n V2 corresponds a vertex of m incident to both fi and /2 , and such a vertex is the origin 
of at least one edge incident to /i with reversal incident to f2- This gives the first equality, and 
the second one follows by symmetry in the indices of the second formula. 

For an element (m, [, t*) G Cg^2, we let f* be the face of m incident to e*, and be the other 
face. Let 

■= '^m^M > 1 < « < degn,(f*) 

be the list of edges of m incident to f* (this notation should not be confused with the notation 
e(i), which denotes edges of m). Besides the process / defined at (fTUj) we will need two other 
auxiliary processes. First, let zq be the first index i such that G f^. Let 

e'l = fmC^io) = a^^H^^'Ht*)) , 

which is incident to and has the same origin as tig, and 

t', = ipi^\e[), l<z<degji). 

We define a process I' as the concatenation 



I' = Z< .../'deg„(T.) 
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Next, let w be defined as the concatenation 



We also consider the function defined for < t < t{w) by T(t) = i if X]o<i<j'"i^^ (e ) — * ^ 

Eo<j<j Vm(e.)' ^(^('^)) = degm(f*) - 1- 

When (m, [1]) = ^'g,2(q, x, D) and {m,l,t*) = Hg,2(m, [1]) for some (q,x, £>) G Qg,2, it is 
straightforward, though a little tedious to check formally, that I' is alternatively described as 
the process interpolating linearly between the values 

I'ii) = liifUe'i)) - l(e'i) , < i < deg„(7j - 1 , (28) 
and /'(degni(/J) = 0, where 

e'l = o-mVm(e*) , 



and where jo is the first index j such that (/?m(e*) G f^,- One checks in this case that 

io-l 

l(e'i) = E^^ (29) 

i=l 

SO that the function X^^Ll^ + ^' reproduces faithfully the labels of half-edges incident to 
Still in this case, the function w reproduces the labels 1 along the face of m-^ incident to e*, 
while T{t) records the label of the edge of m on which the edge of m-^ explored at time [t\ 
projects. 

Lemma 21 Let (q,x,Z)) G Qg,2 md (m,l,t*) = ^g^2i<i,^,D). Assume that mm{'w{t) : ex(t) G 
f^} is attained for a unique value im ofTt, such thatti^ G f*: Also assume that for some R> 0, 



mm{w{t) : ex(t) G f^,T{t) / i^} > min{'u;(t) : ex(t) G f J + R/2 . 



Then 



diam (medqgf;f2(Xq)) 

< 2 max a; (r, diam {i < t' : w'{c'{i)) < inf w' + R/2}) + R + 2 , 

Proof. Recall that Xq is considered as a geodesic metric space, so that the quantity medqg^^'^ 
makes sense. Assume that x G medqg^^'^^l^q)- We first argue that there exists an element of 
medqg^^'^^l^q) ^ distance at most 1 from x. Indeed, if a; ^ ^(q)) then it must stand on 

an edge of the metric graph Xq that links two vertices y, z. Thus any geodesic path from x to xi 
or X2 must pass through y or z, and by symmetry we may assume dq{x, xi) = dq{x, y)+dq{y, xi). 
If we also have dq{x,X2) = dci{x,y) + dq{y,X2) then by subtracting iiq(y,xi) — dq{y,X2) = D 
and y is the wanted vertex. Otherwise, it holds that (iq(a;, X2) = dq(x, z) + (iq(z, X2) so that 

D = dq(x, xi) - dq{x, X2) = dq{y, xi) - dq{z, X2) + (dq(x, y) - d^{x, z)) . 

As edges have length 1, we have dq{x,y) = 1 — d(^{x,z) > 0, so the quantity in brackets is 
in (—1,1), and the others are integers, so we must have dq{x,y) = dq{x,z) = 1/2. But since 
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\dq{z,X2) — (iq(y, 2:2)1 = 1) this implies that (iq(y,xi) — d^[y,X2) has not the same parity as D. 
This contradicts the fact that li — I2 takes values in 2Z established in Sect. 12.21 Thus, we may, 
and will, assume that x E V^(q) up to increasing the diameter by 2. 

We next argue that there is a vertex y G Vi H V2 at dq-distance at most R/2 from x. If x itself 
is not such a vertex, then it is incident only to edges of £'j(q,x, D) for some i G {1, 2}. Assume 
I = 1 by symmetry and note that with the notations of Sect. 12.2] it holds that li(x) = \2{x). Take 
a chain of edges of q starting from x and ending at X2, along which the function 1 decreases, so 
that this chain is geodesic. It is easy to check that the first element e of this chain which belongs 
to E2{(\, x, D) is such that y = e~ is an element of VinV2. Therefore, the chain from x to y can be 
extended as well to a geodesic from x to xi. Thus D = (iq(a;, xi)—dq{x, X2) = dq^{y, xi)—dq{y, X2) 
and y is an element of of med^j2,2(Xq) n V{q). Since x E c{.gxiX2 V ^ geodesic from x 
to xi and from x to X2, it is easily computed that y is another element of q_gxiX2 distance at 
most R/2 from x. 

Therefore, we may, and will, restrict our attention to x E medqg^^'^2'^^i'~'^2; up to increasing 
the diameter by i? + 2. Letting (m, [1]) = ^'g^2(qi x, D), the corresponding x E V{m) satisfies 
dq{x, xi) + dq{x, X2) = 21(x) — l(xi) — l(x2) < dq(2;i, ^2) + R, and together with ([27|) this yields 

l(x) < min 1 + R/2 . 

Now since x is incident to and f^,, we can find a time t with l(x) = w{t), and ex(t) is incident 
to f^, and one sees after a moment's thought that the last inequality can be rewritten as 

l(x) = w{t) < min{u'(s) : tr{s) e fj + • 

These facts imply by hypothesis that T(t) = im- Moreover, by Lemma O if y is another vertex 
such as X, we have dq{x,y) < 2uj{l^^"^, \p — p'\), where ti^{p)~ = x and ej„(p')~ = V with the 
notations of Section l3.3[ Since x and y are incident to /* and we have w^{Cp) = l(x) — l(ei„(0)) 
and w)^(Cp/) = l(y) — l(ej^(0)) which must be both < minw^ + R/2. Thus the result. □ 

By Theorem [21 along the sequence (a„,n > 0), the spaces a' -V4x« under converge to 

[X, d, fj] in distribution, for some limiting [X, d, fi] with law Jy^g^\ Together with lemma \W\ this 
entails that the marked space [Xq, a^^^'^d^, (x, y)] under Qg^^"''' (dq)//q^(x, y) converges to a space 
[X,d,{x,y)] with law y^'^\dX)iJ,'^'^{d{x,y)). On the other hand, the law of ^ g^2{<l, {x,y), D) 
ier the law Qg^l""^ defin 

Q(^/'^)(dq)/.f (d(x,y))- 
(recall p^ ) equals 

I mW-^M r - ■-Mg,2 ((V2 + x(g))"^|P(q, y))rie-/3(^./2+xto))/a d(m, [, t*)) 
LM^ 2 Idim, I, )) .- ^^^ ^ ^^^^^^ ^ ^(^g))-ie-Pin/2+xm/a) 

here, we take the convention |P(q, (x, y))|~^ = when ^(q, (x,y)) is empty. The random 
variable |I'(q, (x,y))| = (dq(x,y) — 1) V is a continuous function of (m, l,t*), since (fST]) . ([28]) 
and pop and the definition of t/; imply 

io-l 

dq{x,y) = 2 min_ ii;(t) — minZ — min/' — l^^ — 2 . (30) 



under the law Q^gl'^'^ defined by the formula 

#P(q,(x,,))(dI)) ^ Qg,2 (F-^|P(q,(x,y))|-ie-/^'^'^/"d(q,(x,y),D)) 
, ™q ^-^--^^ lP(q,(x,y))| Q,,i(l^q-^e-/3Wa) 
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We let d denote the latter quantity with the —2 omitted, which can be associated with any 
element of Cg^2 (not necessarily in 

Lemma 22 The law of ^"■{m,l,t*) under LM^^^^'^^ converges towards the law CLM^^l defined by 

"I 

CLM^'^^ldK [,r)) = 3- 



^Wf^f^ r _ .CLM,,2 {r^h-'e-^^' d(m, l,t*)) 



CLMg,i(r(^ie-/3->) 

Proof. Let H be bounded continuous on Cg 2- Let /ig be on [0,e/2], equal to 1 on [e, oo) and 
interpolate linearly in-between. Let i/e(m, l,t*) = H{m, 1, t*)h£{d)h^{Ti). Then a straightforward 
use of Theorem [5] entails that LI\/1^^2'^^(i?£o^") converges to CLM^'^g (-^^e)- When H is nonnegative, 

this quantity converges as e — > to CLM^^g (if), and in particular CLM^'^g is a sub-probability 
measure. Thus, by dominated convergence, this last convergence holds for H not necessarily 
non- negative. 

On the other hand, we know that (T[/2-|-x(fi'), f ) has same law under LM^^l""^ as {Vq, dq(x, y) + 
2) under Q'f^''\dq)fi'^^{d{x,y)). So {a-Hn/l + x{g)),a-^/^d) under LM^'J'"^ converges in dis- 
tribution to {V,d{x,y)) under J>^g^\d{X,V))fj.'^^{d{x,y)) (here we used Proposition [TO] for the 
convergence of the distance). Since we know that ^ is a.s. diffuse, it holds that d{x,y) > a.s. 
under the latter measure, entailing the tightness of the variables a~^/^D under LM^^g"^ in (0, oo). 

Since a~^Vq = o,~^{t[/2 + xio)) under VM^^l"^ are tight in (0, oo) as well, this entails that for 
some positive r/(e) — > as e ^ 0, 

V.Mft\Hoi-)-Lu[%''\H,oC)\ < ||ii||oor?(e). 

By first taking e small enough and then letting a — > oo, we conclude that V.ufl''\H o ^") ^ 
CLM^'J (if), as wanted. □ 

As in the beginning of Section 17. 2. H we make the assumption that we are working with 
a probability space {Q,J^,P) under which are defined random variables {qa, {xa,ya), Da) with 
respective laws Q^g'^^"'\dq)n®^{d{x,y))\V{q, ix,y))\~'^#'r,{ci,{x,y)){'iD), such that, along a„,n > 1, 

1. the marked spaces [Xa,a~^^^da, {xa,ya)] converge a.s. in PM2 to [X,d,{x,y)] with law 
TTlyP{dX)m2{X,d{x,y)), where vr^ : PM„ x (0,oo) PM„ is the first projection, 

2. a~^^^Da D a.s., the latter having uniform law on {—d{x,y),d{x,y)) given d{x,y), and 

3. letting 'i' g,2{<ia, {xa,ya), Da) = (ma,[lj) and '^g,2{(la, {xa,ya), Da) = (tiTa, [«,**), we have 
a.s. ^"(nxa, la,tl) (m, M*) where the limit has law 

Lemma 23 The following events are of full measure under CLMg^2- 



The quantity mi{w{t) : tY{t) £ f*} attained for a unique time t, which is not a jump 
time ofT{t) and is distinct from and t{w). 

The quantities diam{0 < s < : w^{cf{s)) < iniw^ + 6} for t ^ E decrease to as 5 [Q. 
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Proof. Under the measure CLMg.2 and conditionally on m, r, {w^,t G E{m)), the processes 
are independent Brownian bridges. Since the minimum value attained by a Brownian bridge 
has a diffuse law, and is attained a.s. once by the path at a time which is neither the starting 
nor the ending time of the bridge, the first point is clear. 

Now conditionally on m, r, {w^, t G E{m)), the paths (c^, z^) are independent Brownian snakes 
under CLMg^2- If Sm is the (unique) time at which attains its minimum, then it follows from 
standard properties of Brownian motion that cf (the infimum process of c^) attains the value Sm 
at a unique time tm with probability 1. This means that cl < c^^ < c^, for every t < tm < t' , and 
by continuity the two bounding terms converge towards the middle one as t, t' ^ tm- Combining 
this with the continuity of and the fact that the value iniw^ is attained only at s^, we obtain 
the result. □ 

We can now finish the proof of Theorem [3l By the first point in Lemma [23] and the a.s. 
convergence of ^"(rria, la,ta) to (m, I, t*), we conclude that (ma, la,ta) satisfies the hypotheses of 
Lemma El] for R in the form 6a^^^ with 5 small enough, and every large enough a, almost-surely. 
From the latter lemma and to be able to apply Lemma [TT] to the spaces [Xa,a~^/'^da, {xa,ya)] 
and the delays Da, it suffices to notice that for every rj > 0, there exists 6 £ (0,ry) such that for 
every a, e, 

diam ({z < : wl{cl{i)) < miw^ + a'/^6/2}) < 2arj , 

as follows from the second point of Lemma [23] and the a.s. uniform convergence of the processes 
(j)"" (wl) , {c^) towards w^, c^. Then, taking e > and choosing r] in turn so that 2quj{ip°'{ll),ri) + 
r] < e for every a, e, which is possible by the uniform convergence of V'"(^a) to Lemma [21] 
implies that 

lim sup diam (medqg"~'''''^'"^(a"^/''Xa)) < e. 

Therefore, Lemma [TTI implies that a.s. under 

)(d(X, V))fi^^{d{x, y))dZ^l[_rf(,,j,),,(,,,)] {D)/2d{x, y) , 

all geodesic paths from x to y pass through a unique D-median point. This shows that a.s. 
under the measure ^g^^(d(X, V))/i'^^(d(x, y)), the geodesies from x to y all pass through a 
unique Z)-median point, for Lebesgue-almost every D. Since D-median points form a continuous 
parametrization in D G [—d{x,y),d{x,y)] of geodesies from x to y, this implies that a.s. the 
geodesic between /i-almost all x and y is unique, hence the result. 
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